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We discuss the origin of Yukawa textures in the string-inspired and string derived 
models based on the gauge group SU (4) ® SU (2) l <8> SU (2)r supplemented by a U(l)x 
gauged family symmetry. The gauge symmetries are broken down to those of the 
minimal supersymmetric standard model which is the effective theory below 10 16 GeV. 
The combination of the U(l)x family symmetry and the Pati-Salam gauge group leads 
to a successful and predictive set of Yukawa textures involving two kinds of texture 
zeroes: horizontal and vertical texture zeroes. We discuss both symmetric and non- 
symmetric textures in models of this kind, and in the second case perform a detailed 
numerical fit to the charged fermion mass and mixing data. Two of the Yukawa textures 
allow a low energy fit to the data with a total x 2 of 0.39 and 1.02 respectively, for three 
degrees of freedom. We also make a first attempt at deriving the non - renormalisable 
operators required for the Yukawa textures from string theory. 



1 Introduction 



Over recent years there has been a good deal of activity concerned with understanding 
the pattern of fermion masses and mixing angles within the framework of supersymme- 
try and unification (see next section for a review). The starting point of these analyses 
is the idea that at high energies the Yukawa matrices exhibit a degree of simplicity, 
typically involving texture zeroes, which can be understood as resulting from some 
symmetry. The types of symmetry which have been considered include grand uni- 
fied symmetry to account for the vertical mass splittings within a family, and family 
symmetry to account for the horizontal mass splittings between families. In order to 
restrict the rather ad hoc nature of such models, one may appeal to a rigid theoretical 
structure such as string theory in terms of which the high energy field theory may 
be viewed as an effective low energy supergravity model valid just below the string 
scale. Viewed from this perspective certain classes of unified gauge group and family 
symmetry appear to be more promising than others, and in addition one may hope to 
begin to derive the entries of the Yukawa matrices as low energy non-renormalisable 
operators which arise from the string theory. 

In this paper, guided by the principles outlined in the previous paragraph, we 
investigate the origin of Yukawa textures in a class of models based on the Pati-Salam 
S77(4) x SU{2) L x SU(2) R symmetry with gauged U(l) family symmetries. 

We shall follow both a bottom-up approach, in which the successful textures may be 
extracted from the known quark and lepton masses and quark mixing angles, and a top- 
down approach in which we shall begin to see how the desired operators may emerge 
from a particular superstring construction. This model involves both quark-lepton 
unification, which leads to Clebsch relations to describe the mass relations within a 
particular family, and a U(l)x gauged family symmetry which may account for family 
hierarchies. Thus we are led to vertical and horizontal texture zeroes which are a feature 
of this model. In the earlier parts of the paper we shall focus on something we call 
the string-inspired SU(A) x 0(4) (~ SU(A) <g> SU(2) L ® SU(2) R ) model which contains 
many of the features of a realistic string model such as small group representations and 
a U(l)x family symmetry. Within this simplified model we shall relate the high energy 
textures to the low energy quark and lepton masses and quark mixing angles, and so 
determine by a bottom-up procedure the operators which are likely to be relevant at 
high energies. Later on we shall focus on a particular string construction from which 
we learn how non-renormalisable operators may be generated from first principles. 
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The detailed layout of the paper is as follows: In section ^| we review some ideas 
concerning Yukawa textures, and summarise recent progress in this area. In section ||] 
we briefly review the string-inspired SU(4) x 0(4) model. In section f| we discuss 
symmetric textures in the above model. In section [| we discuss the non-symmetric 
textures. In section |6] we perform a full numerical analysis of the non-symmetric 
models. In section [7] we review the U(l)x family symmetry approach to the model, 
and perform an analysis relevant for the full (symmetric and non- symmetric) model. 
In the subsequent two sections we present a viable string construction of the model 
and indicate how the non-renormalisable operators may arise in the specific string 



construction. Finally section 10 concludes the paper. 



2 Yukawa Textures 



The pattern of quark and lepton masses and quark mixing angles has for a long time 
been a subject of fascination for particle physicists. In terms of the standard model, this 
pattern arises from three by three complex Yukawa matrices (54 real parameters) which 
result in nine real eigenvalues plus four real mixing parameters (13 real quantities) 
which can be measured experimentally. In recent years the quark and lepton masses 
and mixing angles have been measured with increasing precision, and this trend is likely 
to continue in the future as lattice QCD calculations provide increasingly accurate 
estimates and B-factories come on-line. Theoretical progress is less certain, although 
there has been a steady input of theoretical ideas over the years and in recent times 
there is an explosion of activity in the area of supersymmetric unified models. This 
approach presumes that at very high energies close to the unification scale, the Yukawa 
matrices exhibit a degree of simplicity, with simple relations at high energy corrected 
by the effects of renormalisation group (RG) running down to low energy. For example 
the classic prediction that the bottom and tau Yukawa couplings are equal at the 
unification scale can give the correct low energy bottom and tau masses, providing 
that one assumes the RG equations of the minimal supersymmetric standard model 
(MSSM) HQ. In the context of the MSSM it is even possible that the top, bottom and 
tau Yukawa couplings are all approximately equal near the unification scale 0], since 
although this results in the top and bottom Yukawa couplings being roughly the same 
at low energy, one can account for the large top to bottom mass ratio by invoking a 
large value of tan (3 defined as the ratio of vacuum expectation values (VEVs) of the 
1 The next-to-MSSM (NMSSM) with an additional low energy gauge singlet works just as well 
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two Higgs doublets of the MSSM. 

These successes with the third family relations are not immediately generalisable 
to the lighter families. For the remainder of the Yukawa matrices, additional ideas are 
required in order to understand the rest of the spectrum. One such idea is that of 
texture zeroes: the idea that the Yukawa matrices at the unification scale are rather 
sparse; for example the Fritzsch ansatz j4|. Although the Fritzsch texture does not work 
for supersymmetric unified models, there are other textures which do, for example the 
Georgi-Jarlskog (GJ) texture H for the down-type quark and lepton matrices: 
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After diagonalisation this leads to A T = A&, A, 



3A S , A e 



Arf/3 at the scale Mq UT 
which result in (approximately) successful predictions at low energy. Actually the 
factor of 3 in the 22 element above arises from group theory: it is a Clebsch factor 
coming from the choice of Higgs fields coupling to this element. 

It is observed that if we choose the upper two by two block of the GJ texture to be 
symmetric, A i2 = A 2 i, and if we can disregard contributions from the up-type quark 
matrix, then we also have the successful mixing angle prediction 



V u . 



\J Ad/ A s . 



(2) 



The data therefore supports the idea of symmetric matrices, and a texture zero in the 
11 position. Motivated by the desire for maximal predictivity, Ramond, Roberts and 
Ross (RRR) || have made a survey of possible symmetric textures which are both 
consistent with data and involve the maximum number of texture zeroes. Assuming 
GJ relations for the leptons, RRR tabulated five possible solutions for the up-type and 
down-type Yukawa matrices. We list them below for completeness: 
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Solution 3: 



Solution 4: 
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Solution 5: 
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(7) 



where A = 0.22, and the top and bottom Yukawa couplings have been factored out for 
simplicity. These textures are valid at the unification scale. All of the solutions involve 
texture zeroes in the 11 entry. Solutions 1,2, and 4 involve additional texture zeroes 
in the 13=31 positions which are common to both up-type and down-type matrices. 
Solutions 3 and 5 have no texture zeroes which are common to both up-type and down- 
type matrices, apart from the 11 entry. Thus solutions 1,2 and 4 involve rather similar 
up-type and down-type matrices, while solutions 3 and 5 involve very different textures 
for the two matrices. 

Having identified successful texturesQ, the obvious questions are: what is the origin 

of the texture zeroes? and: what is the origin of the hierarchies (powers of the expansion 

parameter A)? A natural answer to both these questions was provided early on by 

Froggatt and Nielsen (FN) [RJ . The basic idea involves a high energy scale M, a family 

symmetry group G, and some new heavy matter of mass M which transforms under 

G. The new heavy matter consists of some Higgs fields which are singlets under the 

vertical gauge symmetry but non-singlets under G. These break the symmetry G by 

developing VEVs V smaller than the high energy scale, There are also some heavy 

fields which exist in vector-like representations of the standard gauge group. The 

vector-like matter couples to ordinary matter (quarks, leptons, Higgs) via the singlet 

Higgs, leading to "spaghetti-like" tree-level diagrams. Below the scale V the spaghetti 

diagrams yield effective non-renormalisable operators which take the form of Yukawa 

2 Over the recent years, there has been an extensive study of fermion mass matrices with zero 
textures Ml. 
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couplings suppressed by powers of A = V/M. In this way the hierarchies in the Yukawa 
matrices may be explained, and the texture zeroes correspond to high powers of A. 

A specific realisation of the FN idea was provided by Ibanez and Ross (IR) 0, 
based on the MSSM extended by a gauged family U(l)x symmetry with 9 and 9 sin- 
glet fields with opposite X charges, plus new heavy Higgs fields in vector representa- 
tions[]. Anomaly cancellation occurs via a Green-Schwarz-Witten (GSW) mechanism, 
and the U(l)x symmetry is broken not far below the string scale ||. By making cer- 
tain symmetric charge assignments, IR showed that the RRR texture solution 2 could 
be approximately reproduced. To be specific, for a certain choice of U(l)x charge 
assignments, IR generated Yukawa matrices of the form: 



3 A \ / z8 z3 z4 \ / z5 z3 



A' 



e e e 

e 3 e 2 e 1 . \ L> 



V e 4 e 1 



e e e 

e 3 e 2 e 1 . A L 



/ e 5 e 3 
e 3 e 
V 1 



\ e 4 e 1 

These are symmetric in the expansion parameters e and e, which are regarded as inde- 
pendent parameters. This provides a neat and predictive framework, however there are 
some open issues. Although the order of the entries is fixed by the expansion param- 
eters, there are additional parameters of order unity multiplying each entry, making 
precise predictions difficult. A way to address the problem of the unknown coefficients 
has been proposed in [IT] where it has been shown that the various coefficients may 



0X1 S6 clS db result of the infra-red fixed-point structure of the theory beyond the Standard 
Model. 

Note that the textures for up-type and down-type matrices are of similar form, 
although the expansion parameters differ. Also note that there are no true texture 
zeroes in the quark sector, merely high powers of the expansion parameter. Thus this 
example most closely resembles RRR solution 2 with approximate texture zeroes in 
the 11 and 13=31 positions. However, without the inclusion of coefficients, the identi- 
fication is not exact. The best fit to RRR solution 2 is obtained for the identification 
e = A 2 , e = A (alternative identifications, like e = A 2 , e = 2 A 3 lead to larger deviations). 
However even this choice does not exactly correspond to RRR solution 2, as can be 
shown by taking solution 2 and inserting the numerical values of the entries: 

1 x KT 4 \ / 5 x 1(T 3 

X u = | 1 x 10" 4 5 x 1(T 2 I , X D = I 5 x KT 3 2 x 1(T 2 2 x 1(T 2 

(9) 



5 x KT 2 1 V 2 x 10" 2 1 



3 The generalisation to include neutrino masses is straightforward [fLOl 
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We compare these numbers to the order of magnitudes predicted by the symmetry 
argument, making the identifications e = A 2 , e = A 
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1 x 10"2 
5 x 10~2 

2 x 1Q- 1 



2 x 10- 3 
2 x 10" 1 
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(10) 

Comparison of Eq.|9] to Eq.|10| shows that while X u is in good agreement, A D differs. In 
Eq.|10|, the 23 = 32 element is an order of magnitude too large. When the unknown 
couplings and phases are inserted the scheme can be made to work. However, some 
tuning of the unknown parameters is implicit. This can be avoided by introducing a 
small parameter 5 into all the elements apart from the 33 renormalisable element, so 
that Eq.|| gets replaced byQ 



(11) 



The idea is that the suppression factor S originates from some flavour independent 
physics, while the parameters e and e control the flavour structure of the matrices. For 
example, suppose we take e = A as in the previous example but scale down the entries 
by a factor of 5 = 0.2. Then we would have, 
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/ 1 x 10- 6 2 x 10- 3 
2 x 10- 3 1 x 10-2 
V 4 x 10~ 4 4 x 10-2 




(12) 



which provides a better description of the numerical values required by the RRR anal- 
ysis for solution 2 in Eq.|^, at the expense of introducing the parameter 5. This ex- 
ample indicates that if family symmetries are to give the correct order of magnitude 
understanding of Yukawa textures without any tuning of parameters, then an extra 
parameter S needs to be introduced as above. 

Another aspect of the fermion mass spectrum that one would like to understand, 

is that of the mass splitting within a particular family. For example the GJ texture 

in Eq.p] provides an understanding of the relationship between the charged lepton and 

down-type quark Yukawa couplings within a given family, and in the simplest U(l)x 

scheme such relations are either absent or accidental, as seen in Eq|| where the form 

4 In our scheme we will have a Unified Yukawa matrix. This, as we are going to see, will imply a 
common expansion parameter for the up and down-type mass matrices and the presence of a factor S 
in the up-quark mass matrix as well. 
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of X E has been fixed by a parameter choice. Unless such parameters are predicted 
by the theory, as in the extension of the initial IR scheme that is discussed in |[L1| , 
the only antidote is extra unification. Then, the leptons share a representation with 
the quarks, and the magic GJ factors of three originate from the fact that the quarks 
have three colours. For example the SO (10) model of Anderson et al [0 (with both 
low energy Higgs doublets unified into a single 10 representation) predicts Yukawa 
unification for the third family, GJ relations for the charged leptons and down-type 
masses, and other Clebsch relations involving up-type quarks. As in the IR approach, 
the approach followed by Anderson et al is based on the FN ideas discussed above. Thus 
for example, only the third family is allowed to receive mass from the renormalisable 
operators in the superpotential. The remaining masses and mixings are generated 
from a minimal set of just three specially chosen non- renormalisable operators whose 
coefficients are suppressed by a set of large scales. The 12=21 operator of Anderson et 
al is suppressed by the ratio (45i/M) 6 , while the 23=32 and operators are suppressed 
by (45b-l/45i) 2 and (45 B _ L S / where the 45's are heavy Higgs representations. In 
a complicated multi-scale model such as this, the hierarchies between different families 
are not understood in terms of a family symmetry such as the the U(l)x of IR. Indeed 
it is difficult to implement a family symmetry in this particular scheme, as the latest 
attempts based on global U{2) |13| show. To be embedded into a string model, GUTs 
such as .SO (10) require k > 1 Kac- Moody levels. With these higher Kac- Moody levels, 
simple orbifold compactifications in which candidate gauge U{l)x family symmetries 
are present do not easily emerge. Nevertheless there has been some progress in this 
direction and three family S'O(IO) and Eq string-derived models have recently been 



classified [14|. Here we restrict our discussion to string constructions based on the 
simpler k = 1 level of Kac-Moody algebras, which are more "string friendly" . 

The SU(4:) <g> 0(4) string model can be viewed as the simplest string- friendly unified 
extension of the standard model which can lead to Clebsch relations of the kind we 



desire. The Pati-Salam gauge group [[15] may be broken without adjoint representations 



and was considered as a unified string model |lg,[[L7j] some time ago. This model has 
recently been the subject of renewed interest from the point of view of fermion masses 
ll8| , and an operator analysis has shown that it is possible to obtain desirable features 
such as Yukawa unification for the third family, and GJ type relations within this 
simpler model. A particular feature of the published scheme which we would like to 
emphasise here is the idea of Clebsch texture zeroes which arise from the group theory 
of the Pati-Salam gauge group. These Clebsch zeroes were used to account for the 
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lightness of the up quark compared to the down quark, for example ]T8|] . However the 
operator analysis of |18|] did not address the question of the hierarchy between families 
(no family symmetry was introduced for example), nor the question of the origin of 
the non-renormalisable operators. Here we shall introduce a U(l)x gauge symmetry 
into the model and combine it with the Clebsch relations previously used, to provide a 
predictive scheme of fermion masses and mixing angles. We shall also ensure that we 
obtain the correct order of magnitude for all the entries of the Yukawa matrices from 
the symmetry breaking parameter, using structures like that of Eq.|TT[ In our case the 
quantity 5 will be identified with a bilinear of heavy Higgs fields which are responsible 
for generating the Clebsch structures, while the parameters such as e will have trivial 
Clebsch structure (singlets under the vertical gauge group) but will generate family 
hierarchies from the flavour symmetry. This corresponds to there being two types of 
heavy Higgs fields: Pati-Salam gauge singlets (corresponding to IR 6 and 6 fields) which 
break the U{l)x family gauge group but leave the Pati-Salam group unbroken, and 
H, H breaking fields whose bilinear forms are U(l)x singlets but transform non-trivially 
under the Pati-Salam gauge group, thereby giving interesting Clebsch structures. The 
non-renormalisable operators of interest must therefore involve both types of Higgs 
fields simultaneously. In view of the unusual nature of such operators, we shall provide 
a string-based discussion of the origin of such operators. 

It is worth emphasising that the main features of the previous analysis, (like the 
assumption of U(l) symmetries, the introduction of singlet fields etc) appear naturally 
in most of the recent string constructions. Therefore, in the final sections of this paper 
we will try to embed our analysis in the context of realistic string models which are 



constructed within the free fermionic formulation |19j of the heterotic string. In doing 
so, we should keep in mind that, in realistic string constructions [^, [21], ^T|] there are 
usually many constraints and in general the resulting field theory is quite complicated. 
Moreover, (in the language of the fermionic strings pl|), within the same choice of 
boundary conditions on the string basis vectors of the world-sheet fermions, there 
are numerous consistent choices of the projection coefficients which result in different 
Yukawa couplings multiplets and the large number of singlet fields which are usually 
present. For this reason we shall try to develop a 'string model' - independent approach 
and begin by considering a field theory S77(4) x 0(4) model, which possesses the salient 
features of a realistic string model and at the same time is simpler to work with. 
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3 The String Inspired 577(4) ® 0(4) Model 



Here we briefly summarise the parts of the model which are relevant for our analy- 
sis. For a more complete discussion see [JIB| . The gauge group is S77(4) (g> 0(4), or 
equivalent ly 

SU(4) ® SU(2) L <g> SU(2) fl . (13) 
The left-handed quarks and leptons are accommodated in the following representations, 

f- = (4,2,l)=(£ £ £ .I)" (14) 

- - ( d R d B d G e + Y 
F- a = (4,l,2) = (^ <, <, 1 J (15) 

where a = 1, . . . , 4 is an SU(4) index, a,x = 1,2 are SU(2) i;R indices, and i = 1, 2, 3 
is a family index. The Higgs fields are contained in the following representations, 

feS = (l,2,2)=(^: (16) 

(where h\ and h 2 are the low energy Higgs superfields associated with the MSSM.) The 
two heavy Higgs representations are 

# Qfe = (4,1,2) =( U i U l f VH ) (17) 
a H a H a H t H 



and 



djj djj d tj c jj 



U H U H U H V H 



H ax = (4, 1, 2) = 

The Higgs fields are assumed to develop VEVs, 

< H >=< v H >~ M GUT , < H >=< P H >~ M G ut (19) 

leading to the symmetry breaking at Mqut 

SU(4) <g> SU(2) L ® SU(2) R — ► SU(3) C <g> SU(2) L ® U(l) y (20) 

in the usual notation. Under the symmetry breaking in Eq.[2(], the bidoublet Higgs field 
h in Eq.|T^ splits into two Higgs doublets hi, h 2 whose neutral components subsequently 
develop weak scale VEVs, 

< h[ >= vi, < h° 2 >= v 2 (21) 
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with tan/3 = vijv\. 



In addition to the Higgs fields in Eqs. [L7])|I1 the model also involves an SU(4) sextet 



field D = (6,1,1) and four singlets 0o and <pi, i = 1,2,3. 0o is going to acquire 
an electroweak VEV in order to realise the electroweak higgs mixing, while tpi will 
participate in an extended 'see-saw' mechanism to obtain light majorana masses for 
the left - handed neutrinos. Under the symmetry property ^1,2,3 — > (—1) X ^1,2,3 and 
H(H) — > (—1) x H(H) the tree level mass terms of the superpotential of the model 
read 0: 

W = X^FiFjh + X 2 HHD + X 3 HHD + X^HF^ + + /jhh (22) 

where /1 =< 0o >~ 0(mw)- The last term generates the higgs mixing between the two 
SM higgs doublets in order to prevent the appearance of a massless electroweak axion. 
Note that we have banned terms which might lead to unacceptably large neutrino- 
higgsino mixing [p2|. The superpotential Eq.|22] leads to the following neutrino mass 



matrix |L6] 



M 



( 1 

rK M GUT ! (23) 




V 

in the basis (vi,i?j,(fik). Diagonalisation of the above gives three light neutrinos with 
masses of the order (m l j) 2 j 'Mqut as required, and leaves right handed majorana masses 
of the order Mqut- Additional terms not included in Eq.^ may be forbidden by 
imposing suitable discrete or continuous symmetries, the details of which need not 
concern us here. The D field carries colour and therefore not develop a VEV but 
the terms in Eq.[22| HHD and HHD combine the colour triplet parts of H, H and 
D into acceptable GUT-scale mass terms |Tj|. When the H fields attain their VEVs 
at Mqut ~ 10 16 GeV, the superpotential of Eq.|22| reduces to that of the MSSM 
augmented by neutrino masses. Below M GUT the part of the superpotential involving 
matter superfields is just 

W = X^QiUjh + \%QiDjhx + XtLiEjhx + A'( /.,//,//, + . . . (24) 

The Yukawa couplings in Eq.|24] satisfy the boundary conditions 

Xi(Mqjjt) = X^Mgut) = X%(M GUT ) = X%(M GUT ) = X%(M GUT ). (25) 



Thus, Eq.(|25|) retains the successful relation m T = at M G jjt- Moreover from the 
relation A^(Mg{/t) = X 1 ^(M G ut), and the fourth term in Eq.^2], we obtain through 
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the see-saw mechanism light neutrino masses ~ 0{rn 2 a / 'Mqut) which satisfy the exper- 
imental limits. 



4 Symmetric Textures 



In this section we briefly review the results of the operator analysis of ref. ||18||) then 
introduce our new approach based on new operators. We discuss the RRR textures as 
a simple example of the new method. 

The boundary conditions listed in Eq.^ lead to unacceptable mass relations for the 
light two families. Also, the large family hierarchy in the Yukawa couplings appears 
to be unnatural since one would naively expect the dimensionless couplings all to be 
of the same order. This leads us to the conclusion that the A 1 / in Eq.|22] may not 
originate from the usual renormalisable tree level dimensionless coupling. We allow 
a renormalisable Yukawa coupling in the 33 term only and generate the rest of the 
effective Yukawa couplings by non-renormalisable operators that are suppressed by 
some higher mass scale. This suppression provides an explanation for the observed 
fermion mass hierarchy. 



In ref. [18] we restricted ourselves to all possible non-renormalisable operators which 



can be constructed from different group theoretical contractions of the fields: 

Oq ~ (Fifyh Ijp) +H.c. (26) 



where we have used the fields H, H in Eqs. |17| , |T8| and M is the large scale M > Mx- 
The idea is that when H, H develop their VEVs, such operators will become effective 
Yukawa couplings of the form hFF with a small coefficient of order Mq UT /M 2 . We 
considered up to n = 2 operators. The motivation for using n = 2 operators is simply 
that such higher dimension operators are generally expected to lead to smaller effective 
couplings more suited to the 12 and 21 Yukawa entries. However, in our field theory 
approach we shall restrict ourselves to the simple case considering only n = 1 operators 
with the required suppression factors originating from a separate flavour sector. We 
will leave the question of the definite origin of the operators for now. Instead we 
merely note that one could construct a FN sector to motivate the operators, or that 
one might expect such operators to come directly out of a string theory. In section |7| 
we shall introduce a U(l)x family symmetry into the model, which is broken at a scale 
Mx > Mqut by the VEVs of the Pati-Salam singlet fields 9 and 9. According to the 
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Table 1: When the Higgs fields develop their VEVs at Mgut, the n = 1 operators utilised 
lead to the effective Yukawa couplings with Clebsch coefficients as shown. We have included 
the relative normalisation for each of the operators. The full set of n = 1 operators and 
Clebsch coefficients is given in Appendix 1. These n = 1 operators were used in the lower 
right hand block of the Yukawa matrices in the analysis of ref . [|l^] . 

ideas discussed in section ^] we shall henceforth consider operators of the form 

(HFf\ ( f) n ft rn \ 
■W){M^) +h - C - (27) 

where M' represents a high scale M' > Mqut which may be identified either with the 
U(l)x breaking scale Mx or with the string scale. We have further assumed the form of 
the operators in Eq.|26| corresponding to n = 1 and glued onto these operators arbitrary 
powers of the singlet fields 9, 9. Note that the single power of {HH) is present in every 
entry of the matrix and plays the role of the factor of 5 in Eq.[ll|. However, unlike 
the previous factor of S, the factor of (HH) here carries important group theoretical 
Clebsch information. In fact Eq.|27] amounts to assuming a sort of factorisation of the 
operators with the family hierarchies being completely controlled by the 9, 9 fields as 
in IR, with m,n being dependent on i,j, and the horizontal splittings being controlled 
by the Clebsch factors in (HH). However this factorisation is not complete since we 
shall assume that the Clebsch factors have a family dependence, i.e. they depend on 
We select the Clebsch factor in each entry from phenomenological arguments. 

As a first example of our new approach we shall consider the RRR textures discussed 
in section 0. Our first observation is that, restricting ourselves to n = 1 operators, there 
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are no large Clebsch ratios between the up-type and down-type quarks for any of the 
operators. This means that it is very difficult to reproduce RRR solutions such as 
solution 2 where the 12 element of the down-type matrix in Eq.|| for example, is 50 
times larger than its up- type counterpart. Of course this can be achieved by requiring 
an accurate cancellation between two operators, but such a tuning of coefficients looks 
ugly and unnatural, and we reject it. On the other hand the n = 1 Clebsch coefficients 
in Table 1 include examples of zero Clebsch coefficients, where the contribution to 
the up-type matrix, for example, is precisely zero. Similarly there are zero Clebsch 
coefficients for the down-type quarks (and charged leptons). The existence of such 
zero Clebsch coefficients enables us to reproduce the RRR texture solutions 3 and 5 
without fine-tuning. Interestingly they are precisely the solutions which are not possible 
to obtain by the standard IR symmetry approach, which favours solutions 1,2 and 4 
and for which the up-type and down-type structures are similar. Thus our approach is 
capable of describing the RRR solutions which are complementary to those described 
by the IR symmetry approach^. To take a specific example let us begin by ignoring the 
flavour dependent singlet fields, and consider the symmetric n = 1 operator texture, 



A 




QM QN 

O w + s.d. N 
N O 



(28) 



where O33 is the renormalisable operator, s.d. stands for a sub-dominant operator 
with a suppression factor compared to the other dominant operator in the same entry. 
Putting in the Clebsch coefficients from Table 1 we arrive at the component Yukawa 
matrices, at the GUT scale, of 



A 



v 



/ 


V 2A^. 


V2\ D 





A 22 



2Af 3 

2-^23 



y/2X D 



12 



12 



xgy/2/VE 




V2X D 



12 



3\&y/2/y/E 




1 

\ 


1 / 



(29) 



(30) 



(31) 



5 In [p3[ , two of us used an alternative approach in order to reproduce the structure of solutions 
1 and 3 of RRR by the implementation of a symmetry. These solutions were found to lead to the 
optimal predictions for neutrino masses and mixings. This has been achieved by a proper choice of 
charges (integer/half-integer) and by imposing residual Z2 symmetries which forbid different entries 
in the up and down-quark mass matrices. 
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where A^ and Af 2 arise from the dominant O^, operator and A^ 2 comes from a sub- 
dominant operator that is relevant because of the texture zero Clebsch in the up sector 
of Oyi- The zeroes in the matrices correspond to those of the RRR solution 5, but 
of course in our case they arise from the Clebsch zeroes rather than from a family 
symmetry reason. The numerical values corresponding to RRR solution 5 with the 
correct phenomenology are, 

/ 5 x KT 3 




a^ = o 3 x lir 3 3 x lir 2 \ , \° = 5 x lir 3 2 x kt 2 j (32) 

V 01 

Thus, the hierarchy A 22 << A^ is explained by a Clebsch zero and a suppression factor 
of the sub-dominant operator. Using Eq.[32] we can read off the values of the couplings 
which roughly correspond to a unified matrix of dominant couplings 

/ 3 x 1(T 3 1 x 1(T 3 \ 
A = 3 x 1(T 3 2 x 10~ 2 2 x 1(T 2 (33) 
V 1 x 1(T 3 2 x 1(T 2 1 / 

where we have extracted the Clebsch factors. We find it particularly elegant that the 
whole quark and lepton spectrum is controlled by a unified Yukawa matrix such as in 
Eq.|33] with all the vertical splittings controlled by Clebsch factors. 

At this stage we could introduce a U{l)x symmetry of the IR kind, and the flavour 
dependent singlet fields in order to account for the horizontal family hierarchy of cou- 
plings in Eq.|33[ In the present case we must remember that there is a small quantity 
S multiplying every non-renormalisable entry as in Eq.|ll], corresponding to the n = 1 
bilinear 5 = ttj which we have required to be present in every non-renormalisable 
entry. Thus we can understand Eq.[33] as resulting from a structure like, 

/ Se 8 6e 3 5e 4 X 

A= Se 3 5e 2 Se (34) 
V 5e 4 5e 1 / 

where we identify e = A = 0.22 and set 5 ~ 0.2 which gives the correct orders of 
magnitude for the entries, rather similar to the case we discussed in Eq.|12| Here 
of course the considerations apply to the unified Yukawa matrix, however, not just 
the down-type quark matrix. The details of the U{l)x family symmetry analysis are 
discussed in section Here we simply note that such an analysis can lead to a structure 
such as the one assumed in Eq ,|3l|. 

A similar analysis could equally well be applied to RRR solution 3. In both cases 
we are led to a pleasing scheme which involves no unnatural tuning of elements, and 
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naturally combines the effect of Clebsch coefficients with that of family symmetry 
suppression, in a simple way. The existence of the Clebsch texture zeroes thus permits 
RRR solutions 3 and 5 which are impossible to obtain otherwise within the general 
framework presented here. 



5 Non- Symmetric Textures 



In this section we up-date the non-symmetric textures based on both n = 1 and 
n = 2 operators introduced in ref . [|T^] , then extend the new approach introduced in the 
previous section to the non-symmetric domain. As in the previous section, we shall 
begin by ignoring the effect of the singlet fields, which will be discussed in section [7]. 

As discussed in Appendix 2 we shall modify the analysis of Ref. [|l^] to only include 
the lower 2 by 2 block Ansatz: 



At 



Of 2 + s.d. 





33 



B R 



Ba 



B 7 



B* 








o l 




QAd 


X 







o 2 




QAd 


X 







o 3 




QAd 


X 







o 1 




QDd 


X 







o 2 




QDd 


X 







o 3 




QDd 


X 







o l 




QMd 


X 







o 2 




QMd 


X 



This is then combined with the upper 2 by 2 blocks considered in ref.[ 

B 2 
B-, 

Ba 



(35) 

(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 
(43) 



where X stands for whatever is left in the 22 position, after the lower 2 by 2 submatrix 
has been diagonalised. The Clebsch coefficients of the n = 2 operators used in EqsjSBj- 
43] are displayed in Table |2] but we refer the reader to ref. |18j for the explicit realisation 
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QUh 2 



1 

V2 

5 







QDhj 
12V2 

1 

5 

V2 



n/5 



9^ 
1 

5 

6 
I 



3^ 
1 

5 







Q Ad 

QDd 
QMd 

o 1 
o 2 
o 3 



Table 2: Clebsch coefficients of n = 2 operators previously utilised. 



of these operators in terms of the component fields for reasons of brevity. The Ansatze 
listed above present problems because of the breakdown of matrix perturbation theory^. 
For purposes of comparison with the new scheme involving only n = 1 operators, we 
will recalculate the predictions for each of the models from ref. [|T^] numerically in the 
next section. 

We now turn our attention to the new approach introduced in the previous section, 
based on n = 1 operators together with singlet fields which for the moment we shall 



ignore. In this case the 21 operator used in ref. |T8| which gave an up Clebsch coefficient 
1/3 times smaller than the down Clebsch is not available if we only use n = 1 operators. 
We must therefore use a combination of two operators in the 21 position that allow 
the up entry to be a bit smaller than the down entry. We require that the combination 
provide a Clebsch relation between A£ and Af x for predictivity. The two operators 
cancel slightly in the up sector, but as shown later this cancellation is ~ 0(1) and 
therefore acceptable. The result of this is that the prediction of V u b is lost; however this 
prediction was almost excluded by experiment anyway, and a more accurate numerical 
estimate which does not rely on matrix perturbation theory confirms that V u b in ref. [18 



is too large. So the loss of the V u b prediction is to be welcomed! The Clebsch effect of 
the 12 operator (with a zero Clebsch for the up-type quarks) can easily be reproduced 
at the n = 1 level by the operator M for example. 

To get some feel for the procedure we will follow, we first discuss a simple example of 
a non-symmetric texture, ignoring complex phases for illustrative purposes. Restricting 
ourselves to n = 1 operators, we consider the lower block to be A% and the upper block 



6 When the magnitudes of H21, H12 and H22 are calculated they are each of the same order in the 
down Yukawa matrix, thus violating the hierarchy in Eq.105 that was assumed in the calculation of 
the predictions. 
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to be the modified texture as discussed in the previous paragraph. Thus we have, 

/ M 
\= M + A O w + s.d. 

V o o c 




(44) 



where O33 is the renormalisable operator. Putting in the Clebsch coefficients from 
Table ^ we arrive at the component Yukawa matrices, at the GUT scale, of 



A 



/ 

A 21 





A 22 



V2Xg/V5 



V2X D 



21 



A 





V2Xg 



Xg/VE 

v/2Af 2 



(45) 



(46) 



(47) 



-3V2Xg/\fE 1 

where A^ 2 and Xg 2 arise from the difference and sum of two operators whose normalisa- 
tion factor of \/5 has been explicitly inserted, and similarly for A 21 and Xg. To obtain 
the numerical values of the entries we use some typical GUT-scale values of Yukawa 
couplings and CKM elements (see ref. |18|) as follows: 



A 



33 



1, A c = 0.002, A s = 0.013, A„ = 0.04, A n 



10" 6 ,A, 



0.0006, A e = 0.0002, (48) 



V cb = 0.05, V u , 



0.22, V, 



ub 



0.004 



(49) 



4.25, tan/5 = 55, m t = WOGeV 



(50) 



where we have assumed, 

a s = 0.115, nib 

The textures in Eqs.^5], ^6] and ^ imply that the 22 eigenvalues are just equal to 
the 22 elements (assuming matrix perturbation theory is valid - see later), and A32 = 
V cb /2 = 0.025. Thus we have = 0.004,Xg = 0.03. The remaining parameters are 



determined from the relations, 



A, 



0, A 



3A e — XgV^X 1 ^/ X s , V u b — X^Vcb/Xc 



(51) 



Note that the up quark mass looks like it is zero, but in practice we would expect some 
higher dimension operator to be present which will give it a small non-zero value. We 
thus have three equations and three unknowns, and solving we find X\ 



2xl0" 4 , \g 
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2 x 1CT 3 , Af^ = 3 x 10 -3 . The difference between A^ and A^ requires suppression of A 
caused by the Clebsch zero in the dominant operator M . Thus the unified Yukawa 
matrix involves operators with the following approximate numerical coefficients, 



3 x 1(T 3 \ 
A = i 3 x 1(T 3 1.5 x 10~ 2 
2.5 x 1(T 2 1 / 



(52) 



where we have extracted the Clebsch factors, and the 22 and 21 values in Eq.|52| refer to 
each of the two operators in this position separately. The numerical values in Eq.|52| are 



not dissimilar from those in Eq.33, in particular the upper 2 by 2 block is symmetrical 
with the same values as before. In this case the lower 2 by 2 block has a texture zero 
in the 23 position, as well as the 31 and 13 positions, but otherwise the numerical 
values are very similar to those previously obtained in Eq.[33} Thus this particular 
non-symmetric texture can be described by a structure of the kind, 

/ 5e bi 3 5e 3 5e bi 9 \ 
A = 5e 3 5e lor2 5e bi ° (53) 
V 5e U9 5e 1 / 

where we identify e = A = 0.22 and set 5 ~ 0.1 as before. Can such a structure for the 
e's be obtained from the U(l)x symmetry? This will be discussed in section |7]. 

There is no reason to restrict ourselves to non-symmetric textures with a zero in 
the 13 and 31 position, as assumed in ref.|l^]. For example the following texture is 
also viable, amounting to a hybrid of the symmetric case considered in Eq.|28] and the 
non-symmetric lower block just considered. 

/ M N \ 
A= M O w (54) 

V o N o c o 33 J 

Here, O33 is the renormalisable operator. We now perform a general operator analysis 
of the non-symmetric case, assuming n = 1 operators for all non-zero entries (apart 
from the 33 renormalisable entry). In this general analysis there are two classes of 
texture: those with universal texture zeroes in the 13 and 31 position (essentially 
n — 1 versions of the textures considered in ref. ||18p and new textures with non-zero 
entries in the 13 and/or 31 position. For now we will not consider the cases with 
operators in the 13 or 31 positions for reasons of brevity. In the general analysis we 
repeat the above procedure, being careful about phases, and obtain some numerical 
estimates of the magnitude of each entry which will be explained in terms of the U(l)x 
family symmetry as discussed in the next section. 
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With the above discussion in mind, we consider the new scheme in which the 



dominant operators in the Yukawa matrix are O33, 0^ 2 , O 



22 5 



021,021 and O12, where 
the last three operators are left general and will be specified later. We are aware from 



A 



(55) 



the analysis in ref.||18|| that O12 must have a zero Clebsch coefficient in the up sector. 
A combination of two operators must then provide a non-zero O21 entry to provide a 
big enough V u b, an additional much more suppressed operator elsewhere in the Yukawa 
matrix gives the up quark a small mass. At Mqut therefore, the Yukawa matrices are 
of the form 

H l2 e i(i,12 x[ 2 

tfaia&e^ 1 + H %1 x{\S^ H 22 x I 22 e i ^ 

H-^x^e^ 32 H 33 e i<f > 33 

where only the dominant operators are listed. The I superscript labels the charge 
sector and xL refers to the Clebsch coefficient relevant to the charge sector / in the 
ij th position. <f>^ are unknown phases and Hij is the magnitude of the effective di- 
mensionless Yukawa coupling in the ij th position. Any subdominant operators that we 
introduce will be denoted below by a prime and it should be borne in mind that these 
will only affect the up matrix. So far, the known Clebsch coefficients are 



X 12 
x 22 



X 



32 






1 



x 22 



X 



D 

32 



= 1 
-1 



r E 

x 22 



X 



32 



-3 
-3. 



(56) 



We have just enough freedom in rotating the phases of ^1,2,3 and ^1,2,3 to get rid of all 
but one of the phases in Eq.[55|. When the subdominant operator is added, the Yukawa 
matrices are 



A 



where we have defined 



ttD,E 
H 2l 











Hge^i H 22 'e i<t>22 ' 

H w x& 

H 12 xf 2 

-^22^22 
H 32 x§2 

H 12 xf 2 
H 22 x 22 
H 32 x 32 






H33 



H 21 




H 21 





H33 




H33 



n 2 ix 21 e 
H 21 x 2 {-< 



I IT Z.U „i(f>21 

+ ii2ix 21 e 
+ -H2i^2i e 



(57) 



(5? 
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We may now remove 22 ' by phase transformations upon Fi,2,3 but 0^ may only be 
removed by a phase redefinition of ^1,2,3, which would alter the prediction of the CKM 
matrix Vqkm- Thus, $\\ is a physical phase, that is it cannot be completely removed 
by phase rotations upon the fields. Once the operators 2 i, 2i , Oi 2 have been chosen, 
the Yukawa matrices at Mqut including the phase in the CKM matrix are therefore 
identified with Hij, H'22, 4>2i- 



6 Numerical Analysis of Masses and Mixing Angles 
from Non-Symmetric Textures 

In this section we discuss the numerical procedure used to analyse the non-symmetric 
cases introduced in the previous section. We shall perform an analysis on the new 
approach based onn = l operators only, and also re-analyse and up-date the original 



scheme of ref. |18[ for comparison. 

The basic idea is to do a global fit of each considered Ansatz to m e ,m M , m u , m c , 
m t , rrid, m s , nib, ®s(M z ), \V u b\, an d \V US \ using m T as a constraint. We use 
the approximation that the whole SUSY spectrum of the MSSM lies at Msusy = frit 
and that the MSSM remains a valid effective theory until the scale Mqut = 10 16 
GeV. Not wishing to include neutrino masses in this analysis, we simply set the right- 
handed Majorana neutrino mass of each family to be 10 16 GeV so that the neutrinos 
are approximately massless and hence their masses do not affect the RGEs below 



M GUT . Recall the parameters introduced in Eq.[57]: <$L = (ft, = H 2 i', = H 2 



J 21 = V, 11 21 = n i\ , -n 2 i = -"21, 
H22 1H22, H u , H 32 ,H 33 . The values of these 8 parameters plus as at the GUT scale 
are determined by the fit. 

The matrices X 1 are diagonalised numerically and \V u i,(Mgut)\, \V us (Mqut)\ are 
determined by 

Vckm = V UL V D l (59) 

where V\jl-,Vdl are the matrices that act upon the (u, c, t)i and (d,s,b)L column 
vectors respectively to transform from the weak eigenstates to the mass eigenstates 
of the quarks. We use the boundary conditions a\(M GUT ) = a 2 (M Gl / T ) = 0.708, 



motivated by previous analyses based on gauge unification in SUSY GUT models [24 



X u ,c,t,d,s,b,e,Li,T, \V US \ and \ V u b\ are then run[| from Mqut to 170 GeV~ m t using the RGEs 

7 A11 renormalisation running in this paper is one loop and in the MS scheme. The relevant renor- 
malisation group equations (RGEs) are listed in ref. . 
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for the MSSM. Below Mqut the effective field theory of the Standard Model allows 
the couplings in the different charge sectors to split and run differently. The Aj are 
then evolved to their empirically derived running masses using 3 loop QCDcgil loop 
QED ||18|| . m e T and A^(m T ) then| fix tan/3 through the relation [[TI| 



sj2m%{m T ) 

cos/3= — , (60) 

vXt{a t ) 

where v = 246.22 GeV is the VEV of the Standard Model Higgs. Predictions of the 
other fermion masses then come from 



vsinp 

v cos (3 



in',..,, ~ A (/ .,./,( /;/', , ,, ; ^_ 



~ ^(WiJ^p (61) 



where m] = 1 GeV. There are twelve data points and nine parameters so we have 
three degrees of freedom (dof). The parameters are all varied until the global x 2 /dof 



is minimised. The data used (with ler errors quoted) is |25 



m e = 0.510999 MeV 

m M = 105.658 MeV 

m T = 1.7771 GeV 

m c = 1.3 ± 0.3 GeV 

m Phys = 180± i2GeV 

m d = 10 ± 5 MeV 

m s = 200 ± 100 MeV 

m b = 4.25 ±0.1 GeV 

= (3.50±0.91)10^ 3 

\V US \ = 0.2215 ±0.0030 

a s (M z ) = 0.117 ±0.005 (62) 

\V c b\ is fixed by if 32 which does not influence the other predictions to a good approx- 
imation and so \V c t,\ and H 32 effectively decouple from the fit. We merely note that 

8 The superscript e upon masses, mixing angles or diagonal Yukawa couplings denotes an empirically 
derived value, whereas the superscript p denotes the prediction of the model for the particular fit 
parameters being tested. 
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in all cases, to predict the measured value of \V c b\, H32 ~ 0.03. Note that no errors 
are quoted upon the lepton masses because m T is used as a constraint on the data 
and because m e ,m M were required to be satisfied to 0.1% by the fit. In this way we 
merely use the lepton masses as 3 constraints, using up 3 dof. We did not perform 
the fit with smaller empirical errors on the lepton masses because of the numerical 
roundoff and minimisation errors associated with high x 2 values generated by them. 
Also, 0.1% is a possible estimate of higher loop radiative corrections involved in the 
predictions. Note that no other theoretical errors were taken into account in the fit. 



The largest ones may occur in derivations of rrib due to the large coupling |26| and 
the non-perturbative effects of QCD near 1 GeV. It is not clear how to estimate these 
errors since the error on m;, depends upon soft parameters which depend on the SUSY 
breaking mechanism in a very model dependent way and non-perturbative QCD is an 
unsolved problem. The correlations between the empirical estimations of the current 
quark masses are also not included. A potentially large error could occur if the ansatze 
considered are not exact but are subject to corrections by higher dimension operators. 
We discuss this point further in section |f[ 

The results obtained from this analysis are given in Table 3. Out of 16 possible 
models that fit the texture required by Eqs. [56|j55| , 11 models fit the data with x 2 /dof< 
3. Out of these 11 models, 5 fit the data with x 2 /dof< 2 and these are displayed in 
Table 3. The operators listed as Ou, O21, O21 describe the structure of the models and 
the entries H22, Hi 2 , H 2 i, cos (ft, H 33 , H22 , H 2 \ are the GUT scale input parameters of 
the best fit values of the model. The estimated la deviation in ag(M^) from the fits is 
±0.003 and the other parameters are constrained to better than 1% apart from cos0, 
whose la fit errors often cover the whole possible range. Out of the predictions shown 
in Table 3, m d discriminates between the models the widest. as(M z ) takes roughly 
central values, apart from model 5 for which the best fit is outside the la errors 
quoted in Eq.|6~2| on as(Mz). m s , \V u b\ are within la of the data point and m c , \V US \ 
are approximately on the central value for all 5 models. Models 3,1 and 4 are very 
satisfactory fits to the data with x 2 /dof< 1. We conclude that the \ 2 test has some 
discriminatory power in this case since if all of the models were equally good, we would 
statistically expect to have 11 models with x 2 /dof< 1, 3 models with x 2 /dof= 1 — 2 
and 2 models with % 2 = 2 — 3 out of the 16 tested. 

We now briefly return to the original models with upper blocks given by Bi_ s 
in Eqs. |36| - P3| |T8| . After again isolating the only physical phase to A^, a numerical 
fit analogous to the above was performed using the same data in Eq.|62]. The main 



22 



Model 


1 


2 


3 


4 


5 


12 ^ 


QM 


QW 


QH 


QH 


QH 


2 i + 2 i 


M + Q A 


G + Q H 


qM + qA 


G + Q H 


R + S 


H 22 /10-' 2 


2.88 


2.64 


2.69 


2.67 


6.15 


Hn/10- 3 


2.81 


4.41 


2.13 


0.70 


1.21 


H 21 /10- 3 


1.30 


5.97 


1.76 


4.33 


1.91 


COS0 


0.87 


1.00 


0.20 


1.00 


0.61 


-H33 


1.18 


1.05 


1.05 


1.07 


4.6 


H 22 '/10- 3 


1.91 


1.87 


1.87 


1.87 


2.87 


H 2l '/10- 3 


1.94 


1.62 


1.63 


1.66 


0.76 


a s (M z ) 


0.119 


0.118 


0.118 


0.118 


0.126 


m d /MeV 


6.25 


1.03 


8.07 


4.14 


11.9 


m s /MeV 


158 


150 


154 


152 


228 


m c /GeV 


1.30 


1.30 


1.30 


1.30 


1.30 


m b /GeV 


4.24 


4.25 


4.25 


4.25 


4.13 


m p t hys /GeV 


182 


180 


180 


180 


192 


\v us \ 


0.2211 


0.2215 


0.2215 


0.2215 


0.2215 


|K 6 |/io- 3 


3.71 


3.51 


3.50 


3.52 


3.50 


tan/5 


59.5 


58.3 


58.3 


58.5 


65.7 


x7 dof 


0.34 


1.16 


0.13 


0.55 


1.84 



Table 3: Results of best-fit analysis on models with n = 1 operators only. Note that 
the input parameters Hij, Hi/, cos 4> shown are evaluated at the scale Mqut- All of the 
mass predictions shown are running masses, apart from the pole mass of the top quark 10 , 
m phys ^ m< (i _|_ 4ag(m t ) ^ rpj^ Qj^jy[ ma t r ix element predictions are at Mz- 
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Model 


Bi 






Ba 


Be, 


Br 


B 7 


Ba 


C(q(M<7) 


0.123 


0.123 


0.123 


0.124 


0.123 


0.124 


0.125 


0.124 


m d /MeV 


7.58 


9.12 


4.64 


6.18 


7.49 


3.63 


3.53 


4.53 


m s /MeV 


215 


240 


179 


210 


217 


179 


200 


187 


mJGeV 


1.29 


1.38 


1.35 


1.16 


1.29 


1.32 


0.86 


1.31 


m b /GeV 


4.19 


4.17 


4.19 


4.19 


4.19 


4.18 


4.20 


4.19 


m p t hys /GeV 


188 


189 


189 


189 


188 


189 


190 


189 


\v us \ 


0.2212 


0.2213 


0.2214 


0.2212 


0.2212 


0.2215 


0.2212 


0.2214 


IK&I/KT 3 


4.52 


4.37 


4.05 


4.22 


4.56 


3.74 


3.85 


3.98 


tan/3 


63.2 


63.6 


63.4 


63.7 


63.2 


63.8 


64.3 


63.6 


X 2 /dof 


0.95 


0.96 


1.00 


1.05 


0.97 


1.16 


1.87 


1.04 



Table 4: Predictions of best-fit analysis on models from ref. [figf ! with n = 2 operators 
included. All of the mass predictions shown are running masses, apart from the pole mass of 
the top quark. The CKM matrix element predictions are at Mz- 



difference in the fit with these models is that there are now 4 degrees of freedom in 
the fit (since there is one less parameter). All eight models in question fit the data 
with x 2 < 2 and these are displayed in Table f|. We do not display the best fit input 
parameters because they are largely irrelevant for the discussion here, la fit deviations 
of as{Mz) are again 0.003 for Bx_ 8 . Note that whereas these models are able to fit 
\V us \,m s ,md,mb,m c fairly well, their predictions of ag{M z ) are high and outside the 
la empirical error bounds. \V u b\ is naturally high in these models (as found in ref. ||18||) 
and this forces as[Mz) to be large, where \V u b\ may decrease somewhat. To fit mj with 
a high as{Mz) requires a large if 33 element and this is roughly speaking why m\ hys 
is predicted to be quite high. In each model the high value of as(M z ) required is the 
dominant source of y 2 apart from B 7 , where m c is low. 

In comparison to the new scheme with n = 1 operators only, the old scheme with 
n = 2 operators fits the data pretty well, although not quite as well as models 1,3,4. 
The old scheme also has one more prediction than the new one. However, the preferred 
models are the ones incorporating the U(l)x symmetry since they go deeper into the 
reasons for the zeroes and hierarchies in the Yukawa matrices. 
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7 U(l) x Family Symmetry in the 577(4) x 0(4) Model 



In our discussion of the symmetric textures, we assumed that we could obtain the same 
structure as IR. Of course, as we have already mentioned, the case we are examining is 
different in two aspects: (a) the fermion mass matrices of the different charge sectors 
have the same origin, and thus the same expansion parameter and (b) all differences 
between these sectors arise from Clebsch factors. As a starting point, we will therefore 
briefly repeat the IR analysis for symmetric mass matrices in our framework; we 
then go on to consider the non-symmetric case, with the goal of being able to reproduce 
the numerical values (at least to an order of magnitude) of the successful ansatze given 
in the previous section. 

The structure of the mass matrices is determined by a family symmetry, U(l)x, 
with the charge assignment of the various states given in Table ||. The need to preserve 





Q % u\ d\ U e\ v\ hi h 2 H H 


U(l)x 


ai ai ai ai ai ai —2a% —2a% x -x 



Table 5: U(l)x charges assuming symmetric textures. 



SU(2) L invariance requires left-handed up and down quarks (leptons) to have the 
same charge. This, plus the additional requirement of symmetric matrices, indicates 
that all quarks (leptons) of the same i-th generation transform with the same charge 
ai. Finally, lepton-quark unification under SU(4) ® SU{2) L <g> SU{2) R indicates that 
quarks and leptons of the same family have the same charge (this is a different feature 
as compared to IR, where quarks and leptons of the two lower generations have different 
charges under the flavour symmetry). The full anomaly free Abelian group involves 
an additional family independent component, U(1)fi, and with this freedom U{l)x is 
made traceless without any loss of generality^. Thus we set a± = — (a 2 + a%). Here 
we consider the simplest case where the combination HH is taken to have zero charge. 
This is consistent with our requirement that it plays no role in the mass hierarchies, 
other than leading to a common factor S for all non-renormalisable entries. 

If the light Higgs h 2 , hi, responsible for the up and down quark masses respectively, 

arise from the same bidoublet h = (1, 2, 2), then they have the same U(l)x charge so 

11 Since we assume that the 33 operator is renormalisable, the relaxation of the tracelessness condi- 
tion does not change the charge matrix since any additional FI charges can always be absorbed into 
the Higgs hi charges. 
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that only the 33 renormalisable Yukawa coupling to h 2 , h\ is allowed, and only the 
33 element of the associated mass matrix will be non-zero. The remaining entries 
are generated when the U(l)x symmetry is broken. This breaking is taken to be 
spontaneous, via Standard Model singlet fields, which can be either chiral or vector 
ones; in the latter case, which is the one studied in IR, two fields 9, 9, with U(l)x 
charge -1, +1 respectively and equal VEVs are introduced. When these fields get a 
VEV, the mass matrix acquires its structure. For example, the 32 - entry in the up 
quark mass matrix appears at 0(e) because U(l) charge conservation only allows the 
term c c th 2 (e/M 2 ) a2 ~ a3 for o 2 > a 3 , or c c th 2 {6 / M 2 ) az ~ a2 } for o 3 > a 2 . Here e = (< 9 > 
/ M 2 y a2 ~ a '^ where M 2 is the unification mass scale which governs the higher dimension 
operators. In IR, a different scale, M 1; is expected for the down quark and lepton mass 
matrices. 

In our case however, all charge and mass matrices have the same structure under the 
U(l)x symmetry, since all known fermions are accommodated in the same multiplets 
of the gauge group. The charge matrix is of the form 

— 2a 2 — 4^3 — 3«3 — a 2 — 2a% \ 

-3o 3 2(o 2 - o 3 ) cx 2 - a 3 (63) 
— a 2 — 2o 3 0:2 — 03 J 

Then, including the common factor 5 which multiplies all non-renormalisable entries, 
the following pattern of masses is obtained (for vector-like singlets): 

/ <5 e |2+6a| £ £ |3a| £ e |l+3a| \ 

V 5el 1+3a l 5 e 1 J 

whereQ a = 03/(02 — 0^3) • We emphasise that the entries in Eq.|>2] describe the mag- 
nitudes of the dominant operators, and do not take the Clebsch zeroes of the different 
charge sectors into account. Note the existence of a single expansion parameter, for 
all three matrices. Another interesting point is that a unique charge combination a 
appears in the exponents of all matrices, as a result of quark-lepton unification. Actu- 
ally, unlike what appears here, in most schemes the lepton mass matrix is described in 
the generic case by two parameters. For a = 1, one generates the structure in Eq. |34| 
for the unified fermion mass matrices. 

Before passing to the non-symmetric case, let us make a few comments on the 

possibility of having chiral or vector singlets, as well as on the charge of the Higgs 

12 In this simplest (and more predictive) realisation, hi, « h t therefore we are in the large tan f3 
regime of the parameter space of the MSSM. 
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fields. Suppose first that 9 is a chiral field. From the form of the charge matrix, we 
observe that if the 22 and 23 entries have a positive charge, a 3 is negative (for all these 
entries to be non- vanishing at the same time). Moreover the hierarchy 1:3 between the 
23 and 12 elements indicates that a 2 would have to be zero in the chiral case, and thus 
the 13 element would tend to be larger than desired. We can say therefore that in the 
symmetric case with vector fields generates the mass hierarchies in a more natural way. 

Concerning the h%, h 2 higgses, there are two kinds originating from free fermionic 
string models: those coming from Neveu-Schwarz sector which in general have integer 
(including zero) U{l)x charges, and those arising from twisted sectors, which usually 
carry fractional U(l)x charges. Which of these cases acquire VEVs, is decided from 
the phenomenological analysis. For example, to obtain the structure of Eq.|34| we see 
that the charges of ma y n °t be zero, since in such a case the 12 element which is 
proportional to the Higgs charge would be unacceptably large. For the non-symmetric 
case of course this feature does not necessarily hold. Finally, the H, H fields (the SU(4) 
higgses) tend to be non-singlets under extra U(l)x symmetries. We now proceed to 
discuss the non-symmetric case, which in the framework of U(l)x symmetries has 



been extensively studied in [27fl. Here, we will examine what constraints one may put 



on the various possibilities for non-symmetric textures, in the model under study. 
The charge assignment for this case appears in Table M. Fields that belong to the 



Qi u\ d\ U e\ vt h x h 2 H H 



U(l) x A aci on a>i a { -(3 3 -a 3 -f3 3 -a 3 x -x 



Table 6: U(l)x charges for non-symmetric textures. 

same representation of S77(4) <g) SU(2)l <%> SU(2)r are taken to have the same charge. 

Again, it is clear that all fermion mass matrices will have the same structure. With 

this charge assignment we may proceed as in the symmetric case, and calculate the 

possible mass matrices that may arise. The charge matrix is now 

/ _ ft2 _ 2a 3 -fc- 2(3 3 a 2 -a 3 -(3 2 - 2(3 3 -fa - 2(3 3 \ 

-a 2 - 2a 3 + (3 2 - f3 3 a 2 -a 3 + (3 2 - (3 3 f3 2 - f3 3 (65) 
\ — a 2 — 2a 3 a 2 — a 3 / 

We now want to find which charge assignments may generate a mass matrix as close 

as possible to the form in Eq.^3], keeping in mind that there is no reason to restrict 

ourselves to non-symmetric textures with a zero in the 13 and 31 position. 

In what follows, we will check whether it is possible to generate the hierarchies in 
the effective low energy Yukawa couplings required by our ansatze and the data. The 
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required couplings are detailed in Table 3. Initially, we determine if we can obtain the 
correct structure by chiral singlet fields. We assume for a starting point that for the 32 
entry we have a 2 — a% > (without a loss of generality since we can always choose the 
sign of one entry in the charge matrix). The 23 entry has to be small (it is assumed to 
be zero in the ansatze in Eq.^), indicating that (a) either /3 2 — /3 3 < or (b) f3 2 — P3 
is positive and large (> 2). Case (b) is excluded, since it would indicate that the 22 
charge, which is always the sum of the 23 and 32 charges, would be unacceptably large 
as well (which implies that H 22 < H 32 , in contradiction to the fits in Table 3). What 
about )? A negative number must not dominate the 22 entry in the chiral case, 

thus \j3 2 — /?3 1 would have to be smaller than \a 2 — a&l. This clearly contradicts the 
required hierarchy between the 22 and 32 elements and so the required couplings can 
not be naturally described by a model with only a chiral U(l)x Higgs 9. 

For this reason we are going to look for solutions in the case of vector singlets, where 
it is the absolute value of the charges that matters. Here, the important difference from 
the previous case is that a solution with a small and positive a 2 — a 3 and a large negative 
02 — @3 is allowed. The 23 and 32 elements have the correct hierarchy, while the 22 
element can also be sufficiently small, as a result of a cancellation between terms of 
opposite sign, with the negative contribution being dominant. What can we say about 
the rest of the structure and how restrictive should we be when looking for solutions? 
We could allow for a small asymmetry between the 12 and 21 entries. Actually, Xf 2 
can be slightly larger than A^. This, combined with the fact that there are unknown 
coefficients of order unity indicates that we can have an asymmetry of order e between 
the 12 and 21 entries. We will discuss solutions with such an asymmetry, even in the 
case that \® 2 < X 21 , due to this coefficient ambiguity as well as the ambiguity in the 
experimental value of the up and down quark. We also need not drop solutions with a 
large 13 or 31 entry, if they are compatible with the numerics. 

On this basis, we have looked for solutions in the following way: for the charges 
of the elements 12-21-22-32 we made all possible charge assignments (such that lead 
to a maximum 4 th power in terms of the expansion parameter for the resulting mass 
matrices, for the 12 and 21 entries). This fixes all charges a 2 , /3 2 , $3 each time. 
We then looked at what the charges of the other entries are and whether the generated 
hierarchies are consistent with the phenomenology. 

The restrictions we require in order to identify a viable solution, are (besides of 
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course that the only renormalisable term is in the 33 position) 



charge(ll) 


> 


charge(12) 


charge(ll) 


> 


charge(21) 


charge(21) 


> 


charge(22) 


charge(12) 


> 


charge(22) 


charge(13) 


> 


charge(22) 


charge(31) 


> 


charge(22) 


charge(32) 


< 


|charge(22)|0(e) 


charge(12) 




|charge(21)|0(e) 


charge(23) 


> 


charge(22) 



Then, we end up with the following possibilities: 
Case 1: 

a 2 = -2/3, a 3 = -5/3, fa = -2, fa = 0, Y uM 
Case 2: 

a 2 = -1, a 3 = -2, /3 2 = -2, /3 3 = 0, Y u ,d,t ■ 

Case 3: 

a, = -4/3, a 3 = -7/3, ft = -2, /3 3 = 0, F u> ^ 
Case 4: 

a, = -4/3, a 3 = -l/3, ft = 0, ft = -2, 3^ 
Case 5: 

a 2 = -4/3, a 3 = -7/3, ft = -3, fa = 0, 



5e 6 


5e 3 


5e 2 


5e 2 


5e 


5e 2 


5e 4 


5e 


1 


5e 7 


5e 3 


5e 2 


5e 3 


5e 


5e 2 


5e 5 


5e 


1 


5e 8 


5e 3 


5e 2 


5e 4 


5e 


be 2 


5e 6 


5e 


1 


5e 6 


5e 3 


5e 4 


5e 4 


5e 


5e 2 


be 2 


5e 


1 


5e 9 


5e A 


5e 3 


5e 3 


5e 2 


5e 3 


5e 6 


5e 


1 
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/ 5e 8 Se 4 Se 3 
Se 3 Se 5e 2 
\ 5e 5 5e 1 



/ 5e 10 be 4 Se 3 



(72) 



5e 4 6e 2 5e 3 
V 5e 7 5e 1 



(73) 



Case 6: 

a 2 = -1, a 3 = -2, ft = -7/3, ft = -1/3, F W)d) * 
Case 7: 

a 2 = -5/3, a 3 = -8/3, ft = -3, ft = 0, F U)d) ^ = 
Case 8: 

a 2 = -4/3, « 3 = -7/3, ft = -7/3, #, = -1/3, 
Case 9: 

« 2 = -4/3, a 3 = -l/3, /3 2 = -1/3, ft = -7/3, 



Let us also list for completeness a few cases with a larger splitting between the 21 and 
12 entries (up to 0(e 2 ) ): 
Case 10: 



/ 5e 9 be 4 5e 3 

Se 4 5e be 2 ) (74) 
V Se 6 5e 1 



5e 7 


5e 4 5e 5 


Se 4 


5e 5e 2 


5e 2 


5e 1 



a 2 = -4/3, a 3 = -l/3, /? 2 = 1/3, ft = -5/3, Y uAe = 
Case 11: 

« 2 = -2/3, « 3 = -5/3, ft = -7/3, ft = -1/3, K,^ 



5e 5 5e 2 <5e 3 
5e 4 5e 5e 2 | (76) 
5e 2 5e 1 



/ be 7 5e 4 5e 3 

5e 2 Se 5e 2 | (77) 
V 5e 4 Se 1 



Of course, here we also have the cases with the opposite charge assignment^. 
Among the various choices, we see that 

• The charge of the Higgs fields h± t 2 is always different from zero. 

• There are cases where the 13 and 31 elements are large. 



13 The presence of fractional charges implies the existence of residual discrete symmetries after the 
breaking of the abelian symmetry. 
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We may now examine the results of Table 3 in the context of the XJ(l)x symmetry 
discussion above. We take all models that fit the data with x 2 /dof< 1, i.e. models 
1,3,4. We define in each of these models, H^ p as being the dimensionless and dominant 
effective coupling constants in the SU(4)®SU(2)i®SU(2) fi unified Yukawa matrix for 
the best fit parameters. 



Then, model 1 has 



H, 



emp 
ij 





0.001 




0.003 
0.03 
0.03 1 



(7* 



We see that case 1 above does not fit this pattern very well if all dimensionless couplings 
are ~ 0(1) because in case 1, H 2 ± is suppressed in comparison to H 12 . Cases 4,9 do 
not possess approximate texture zeroes in the 31 position and this would affect \V u t,\ 
strongly. Similar objections can be raised about other cases, except for cases 2,7,8. 
Case 2 with e = 0.21, 5 = .14 yields 



2.10~ 6 

0.001 

6.10~ 5 



0.001 6.10~ 3 
0.03 6.10~ 3 
0.03 1 



(79) 



which fits Eq.[T^ well apart from a factor ~3 in the 12 position. The next sub-dominant 
operator in the 22 position needs to be 2.10 -3 according to Table 3. The values of e and 
5 used in Eq. [T^ give the subdominant operator in the 22 position to be ~ 6.10 -3 . This 
is acceptable, but a closer match occurs for the next higher dimension operator, which 
has magnitude ~ 10~ 3 . An ambiguity occurs in that we have not set the normalisation 
of the sub-dominant operator due to its numerous possibilities and so the original 
discrepancy factor of ~ 3 could easily be explained. Below, we do not consider the 
numerical size of the sub-dominant operator because it is clear that some operator can 
be chosen that will fit the required number well. If the charge assignments under the 
U(l)x symmetry were the same as in this case, we would have succeeded in explaining 
why the assumption of texture zeros was valid. For example, the 13 element in Eq.|79| 
being 6 x 10~ 3 instead of zero only affects mixing angle and mass predictions by a 
small amount. We have also explained the hierarchies between the elements in terms 
of the different mass scales involved in the non-renormalisable operators by not having 
to choose dimensionless parameters of less than 1/3 (or greater than 3). Case 7 with 
e = 0.36, 5 = 0.08 gives 



2.10~ 6 

0.001 

6.10~ 5 



0.001 4.10~ 3 
0.01 4.10~ 3 
0.03 1 



(80) 
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We should note that at this level, we may naively expect 8% corrections to the con- 
straint in Eq.[8(] through the next order of 5 operators in each element. We could have 
attempted to include these possible errors in the numerical fits but we did not due to 
the fact that they are very model dependent. Deeper model building in terms of con- 
structing the non-renormalisable operators out of extra fields or examining underlying 
string models would be required to explain why this should not be the case. It should 
also be borne in mind that explanations for exact texture zeroes can be made in this 
context by setting fractional U(l)x - charges on the heavy fields in the operators, or 
by leaving certain heavy fields out of the FN model. Case 8 with e = 0.36, 5 = 0.08 
gives the same results as in Eqj80|, except with the (22) element as 0.03. 

From Table 3 we see that model 3 (the model that fits the data the best) has 



h: 



emp 





0.002 




0.002 
0.03 
0.03 1 



Choosing e = 0.26, 5 = 0.12 in case 2 gives a good match to Eq.RT| 



9.10- 8 
0.002 
6.10~ 5 



0.002 8.10~ 3 
0.03 8.10~ 3 
0.03 1 



(82) 



Case 7 with e 
match as well. 



0.40, 5 = 0.07 or case 8 with the same e and 5 both give a fairly good 



Model 4 is different in the sense that it possesses a hierarchy between the 12 and 
21 entries of the effective Yukawa couplings: 



H 



emp 



0.0007 
0.004 0.03 
0.03 1 



(83) 



Here, case 1 with 5 = 0.2, e = .15 predicts 



3.10~ 7 0.0007 4.10- 3 
0.004 0.03 4.10- 3 
10- 5 0.03 1 



»4) 



an extremely good match to Eq| 
match also. 



Case 6 with e = 0.28,5 = 0.11 provides a good 



Thus we see that we can explain the hierarchies and texture zero structures of the 
models that fit the data best. In general, it seems likely that we have enough freedom 
in setting charges to attain the required hierarchies for the Yukawa matrices. 
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8 The String Model 



In the following, we will present a semi-realistic string model which provides an exis- 
tence proof of how previously described non-renormalisable operators may be generated 
from first principles using string theory. Before this, let us briefly comment on how 
the basic features of the U(l)x symmetries that we have discussed arise in string con- 
structions. 

In realistic free fermionic string models [19, [L7| there are some general features: 
At a scale M string ~ 5g str ing x 10 17 GeV, one obtains an effective N — 1 supergravity 
model with a gauge symmetry structure which is usually a product of non-Abelian 
groups times several U(l) factors. The non-Abelian symmetry contains an observable 
and a hidden sector. The massless superfields accommodating the higgs and known 
chiral fields transform non-trivially under the observable part and usually carry non- 
zero charges under the surplus £/(l)-factors. The latter, act as family symmetries in 
the way described above. Some of them are anomalous, but it turns out that one 
can usually define new linear U(l) combinations where all but one are anomaly-free. 
The anomalous U(l) is broken by the Dine Seiberg Witten mechanism p8| |, in which a 
potentially large Fayet-Iliopoulos D-term is generated by the VEV of the dilaton field. 
A D-term however breaks supersymmetry and destabilizes the string vacuum, unless 
there is a direction in the scalar potential which is D-flat and F-flat with respect to 
the non-anomalous gauge symmetries. If such a direction exists, some of the singlet 
fields will acquire a VEV, canceling the anomalous D term, so that supersymmetry 
is restored. Since the fields corresponding to such a flat direction typically also carry 
charges for the non-anomalous D-terms, they break all U(l) symmetries spontaneously. 
For the string model in ref . ]17|] , the expected order of magnitude for the VEV of the 
singlet fields is < >~ (0.1 — 0.3) x M string . Thus, their magnitude is of the right 
order to produce the required mass entries in the mass matrices via non-renormalisable 
operators. 

As an application of the above procedure, we will make a first attempt to derive the 
relevant operators for the mass matrices of the model based on the work in ref. [|HJ . The 
string model is defined in terms of nine basis vectors {5, b%, b 2 , b 3 , 6 4 , b 5 , 6 6 , a, (} and a 
suitable choice of the GSO projection coefficient matrix. The resulting gauge group has 
a Pati-Salam (577(4) x SU(2)l x SU(2)r) non-Abelian observable part, accompanied 
by four £7(1) Abelian factors and a hidden 577(8) x U(l) symmetry. 
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In the following, for convenience, we denote a set of complex right fermions with the 
letters j 1 ^ 1 '" 5 , y? 1 '" 6 , fj 123 , z 12 } and real right fermions with {y 1 '" 6 , a) 1 '" 6 }. Now, a specific 
model is defined in terms of a set of boundary conditions on the phases picked up by 
the fermions when parallel transported around non-contractible loops. The model is 
derived from the following basis p 



s = {^,x 12 - 6 ;0-.-o} 

h = {^,X 12 ,Z/ 345 ¥ 456 :* : "V[ 

b 2 = {^, x 34 ,y 1 V 1 W 6 ^^V} 

b 3 = {^, X 56 ,^ 123 V 234 ;^ 1 - 5 r ? 3 } 
b A = {^, X 12 ,y 3 V 6 ,^ 45 

h = {r,x 3 WV"> 15 ^ 15 

h = {0,0, 2/V> 15 u> 15 ;* 1 'YV 2 V} 
C = {0....0 ;^V- 6 } 
a = {y 46 y 46 ,^ 46 ;^ 123 r / 12 z 12 } 

All world sheet fermions appearing in the basis are assumed to have periodic boundary 
conditions, while those not appearing are anti-periodic. An immediate consequence of 
using only periodic and anti-periodic boundary conditions is that the resulting gauge 
symmetry is in general a product of SO(n) groups. Thus, in the above basis for 
example, the complex world sheet fermions tj/ 1 '— ' 5 define an SO (10) symmetry which 
is broken by the last vector a into 50(6) <S> 0(4). Now, bearing in mind that this 
part will be interpreted as the observable gauge symmetry, we observe the isomorphics 
50(6) ~ SU(4), 0(4) ~ SU(2) ® SU(2). The two SU(2)s are going to accommodate 
the left and right components of the matter fields. Thus, the resulting gauge symmetry 
is isomorphic to the Pati-Salam gauge group. Thus the complete symmetry of the model 
under the above choice is 

[SU(A) x SU(2) x SU(2) x U(l) 3 } x [SU(8) x U(l)] h (85) 

where the subscripts (o, h) denote the observable and the hidden part respectively. 



With the specific choice of the projection coefficient matrix in [17|], one obtains three 
chiral families in the (4, 2, 1) + (4, 1, 2) representations of the PS symmetry, and two 
higgs pairs transforming as (4,1,2) + (4,1,2), all arising from the sectors 61,2,3 and 
b 4 , b 5 . 
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In particular, the massless spectrum contains three (^1,3,4)1, = (4,2,1) represen- 
tations obtained from the sectors 61,3,4, which accommodate the left handed fermion 
fields. There are five (4,1,2) representations (^1,4,5, F2, F^)r named after the corre- 
sponding sectors and two H i 5 = (4, 1, 2) arising from the sectors 6 4 , 5 . Thus, two linear 
combinations of the F 4 will play the role of the GUT higgs H, while the remaining 
three F's accommodate the right- handed fermions. The spectrum includes also bidou- 
blets hi = (1,2,2)*, sextets Di = (6, 1,1)* and a sufficient number of singlet fields 
&ij,£j,Ck- A certain number of singlets should develop VEVs in order to satisfy the 
flatness conditions and give masses to unwanted colour triplets and exotic states. 

In addition, one obtains fractionally charged states which arise in non-standard 
representations of the PS- symmetry, namely (1, 1, 2), (1, 2, 1) and one pair (4, 1, 1) + 
(4,1,1). Finally, under the hidden gauge group, one obtains 10 irreducible represen- 
tations Zi, Zi sitting in the 8 of SU(8) while carrying quantum numbers under all five 
U(l) symmetries of the model. All states are divided to those arising from the Neveu- 
Schwarz (NS) and Ramond (R) sectors. In particular the NS-sector gives the graviton 
multiplet as well as the singlet fields $y, sextets and the bidoublets h 3 , h 3 . 



9 Calculation of Tree Level and Non Renormalis- 
able Operators in the String Model 



To calculate the superpotential of the model, one needs to obtain vertex operators 
for all physical states of the theory. To construct vertex operators for the states of a 
given model, every world- sheet fermion has to be represented by a conformal field. 
In the case that a representation of the model can be fully factorized in a left and a 
right moving piece, one can pair them up to bosonized fields. Now, according to the 
definition of the supersymmetry generator S in the above basis of our model, one can 
conclude that the left moving fields \ % can be bosonized (x 1 ± iX 2 )/V% = exp{±z5i 2 } 
and similarly for the x 3,4 an d X 5 ' 6 pairs. N = 1 supersymmetry implies the existence 
of an extra current, which is expressed in terms of S*j as follows |37| 



and which is extended to three f/(l)'s generated by Su, S34, S^q. 

The Yukawa couplings in four dimensional superstring models correspond to expec- 



J(q) = id q (Si2 + S u + S 56 ) 
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tation values of the form 



< J d\ J d% J eZW(9i)^(<Z2)V 3 B (?3) > (87) 



where the Vf' B are the vertex operators for the fermionic (F) and the bosonic (B) 
states, while (71,2,3 are the two dimensional coordinates. Thus, a vertex operator 
for any physical state is a collection of conformal fields that represent the quantum 
numbers of the state under all symmetries of the model. The piece of the vertex 
operators involving the bosonized left moving fields \ l is given for the bosons by 
V B X ~ exp{aSi 2 } exp{/3S 3i } exp{"fS 56 }. Similarly, for the fermions, V F x i 2 ~ exp{( 



a 



1/2)5*12} exp{(/3— 1/2)5*34} exp{(7 — 1/2)556}. The subscripts —1, —1/2 refer to the cor- 
responding ghost numbers. The total ghost number should add up to —2, thus in trilin- 
ear terms the non- vanishing couplings are proportional to the correlator < V F V F V B >. 
In non-renormalizable contributions, the remaining vertex operators Vj 3 ■ ■ ■ V B have to 



be 'picture - changed' in the zero picture [3/]. In general, a particular correlator is non 



- vanishing, only if it is invariant under the three Z7(l)'s. In addition it has to respect 
the usual (right moving) gauge invariance and other global symmetries. For example, 
pure NS- couplings are possible only at tree level. The same is true for higher order 
couplings involving only Ramond fields, and so on. A complete list of rules is found in 



If we imply the well defined set of rules to calculate the Yukawa interactions in the 
present string model, we obtain the following tree-level terms that are relevant to our 
discussion 



+ Cihhi2 + Uh 3 h 12 + $12/112^12 + $3^12^12 

+ 6^12 + 6lMl2 + $12^12^12 + $3^12^12 + ■ ■ • (88) 

where the {• • •} stand for terms involving exotic and hidden fields and other couplings 
irrelevant for our purpose. The F-flatness conditions are derived for the complete tree- 



level superpotential, which is given in |T7j and involves in total 18 singlet fields. Five of 
these fields, namely $i,...,5 have zero quantum numbers under the U(l) groups, while 
the rest of the fields (denoted by £1,2,3,4, £1,2,3,4, Ci,2, Ci,2, $12, $12, $u,) have 
non-trivial quantum numbers. 
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From the above, it is clear that only a few Yukawa couplings are available for fermion 
mass generation at the tree level. The missing terms are expected to be obtained from 
non-renormalisable (NR) terms. In the case of the PS - symmetry we expect NR-terms 
of the form 

FFh etc (89) 

string 

which act as effective mass operators once the fields H, H, and $j j get VEVs. The scale 
where the higgs fields H, H obtain their VEVs is determined from phenomenological 
requirements and renormalisation group analysis|38[] of the particular model. Moreover, 



the singlet VEVs are not completely arbitrary since they should satisfy the D- and F- 
flatness conditions. In general, the D-flatness conditions read 

Y^Qx \< *i >| 2 +^^Tr a ce{Q u{1)x }M 2 Pl = (90) 

£QjJ<$ ?; >| 2 =0 (91) 



where < $j > are the singlet VEVs and g stands for the unified gauge coupling at 
M string . U(l)x in (p0| ) is the anomalous U(l) combination and Q x the corresponding 
U(l)x charge of the singlet $j. Eq. ([H]) holds for all the non- anomalous U(l) 
symmetries of the particular model. From the relations ( |9"0"| , |9~TD , it is clear that the 
order of magnitude of the VEVs of the singlet fields is determined by the Trace - term. 
Thus, we expect that 

< * >2 ~ ° (^m^) M " (92) 

In particular, for the string model in ref. [17]], Tr[Qx] = 72, therefore the order of 
magnitude for the singlet fields is < >~ (0.1 — 0.3) x M string . (See also appendix 3 
for the details.) This indicates that the singlet VEVs have the correct magnitude, in 
order to produce the required mass entries in the mass matrices via the non - renor- 
malisable operators of Eq.flgSp. We also note here that the spontaneous breaking of the 
anomalous U(l) symmetry introduces one more mass scale Mx in the theory, which 
is characterised by the magnitude of the related singlet VEVs. Thus, one naturally 
expects the hierarchy M str i ng > Mx > Mqut- 

One possible choice of non-zero VEVs is 

< >,< $12 >,< 6 >,< 6 > (93) 

and < Z 5 >,< Z' 8 >^ of the hidden fields. Solving the flatness conditions (ap- 
pendix 3), one finds that the order of magnitude of the singlet VEVs is J 3 ^ in Planck 
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units. It is easy to see that the choice (|93|) satisfies trivially the F flatness condi- 
tions. We should point out however that this choice is not unique. There are other 
cases which also satisfy conditions (|90|j9l|) and hopefully, a solution which meets the 
phenomenological requirements does exist. The non- zero VEVs in (^) provide all 
dangerous colour triplets with masses from tree-level superpotential terms. Here we 
would like to investigate if they are also capable of producing the relevant operators for 
the fermion masses. This computation will prove to be a rather hard task mainly due 
to the rapidly increasing number of NR-operators as the calculation proceeds to higher 
orders. We will see however, that the pattern of the fermion mass matrices described 
in the previous sections is basically obtained. 

We will first start the examination of the tree level superpotential. Due to the 
string symmetries and the U(l) charges of the superfields, as can be seen from (j88"|) 
only three terms relevant to the fermion masses exist at three level 

1 



W — > F^F^hu + — ^F^rF^r^ + F 3 nF 3 Lh 3 



(94) 



Here, hi 2 ,h 3 are bidoublets, (2 is a singlet, while the Fl,r chiral fields have been 
presented previously. We may give a non zero VEV to one of the two bidoublet higgs 
fields (or to a linear combination cos#/i 12 + sm6h 3 ) and support one generation with 
masses at the tree level. Since there are more than one doublets in the spectrum, first, 
we should determine the massless state along the chosen flat direction. At the tree 
level, the bidoublet higgs mass matrix obtained from the relevant terms is 



(h 3 ,h 12 ,h 3 ,h 
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(95) 



with ip = $3/2. In order to have at least one non-zero eigenvalue we impose det[m^] = 
(£1^4 — £i£/0 2 = 0, which is satisfied for any value of the $3, $12, <&u VEVs, provided 
£1^4 = ^1^4. The choice ([JB|) is consistent with these requirements. Moreover, it leaves 
h 3 , h%2 massless at three level. We then let h\2 develop a VEV and give masses to the 
top, bottom and tau particles living in the F^l, F^r representations. The h 3 bidoublet 
is expected to receive a mass from a NR- term. Thus, to proceed further, we need 
the contributions of the non-renormalisable terms. As in the tree level case, a non- 



vanishing NR- term of the superpotential must obey all the string selection rules [37| 
and be invariant under all the gauge and global symmetries. Since here we discuss 
the fermion masses, we are primarily interested in those operators contributing to the 
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corresponding matrices. At fourth order, we find no relevant terms. At fifth order, 
there are several operators which in principle could contribute to the fermion mass 
matrices. We list them here, 

FsrFM^ , F 3R F 3L h 3 <S>l A5 (96) 
F 1R F 1L h 12 £ 2 <S> 2 , F 5R F 4L h 12 $l 2 (97) 
F5 R F4LF 2R F 2R hi 2 (98) 

(scaled with the proper powers of M string ). Let us analyse the above contributions in 
terms of the particular flat direction chosen here. It is clear that, irrespective of the 
choice of the singlet VEVs, the terms (|96|) do not add a new contribution since they 
constitute small corrections to the already existing tree-level term F 3 F 3 h 3 . Moreover, 
within the given choice of our flat direction, < $j >=< ( 2 >= 0, the terms ( [97] ) do 
not also generate any new fermion mass term. Thus, there is only one term which 
contributes to the fermion mass matrices, namely the operator of Eq. (|98"D . This is 
a n — 1 operator according to our classification in the previous sections. We have 
already interpreted F 5R and F^ as the left and right components of the third fermion 
generation. Up to now we have not determined which of F 2R , F^ is going to play 
the role of the second family. The 5 th order operator still leaves this undetermined 
since both fields enter in the operator in a symmetric way. Thus there are two options. 
Either we set < F 2R >= 0, or we have to rename the fields so that F 3L ,F 3R are the 
3 rd generation fermions and h 3 the massless higgs. In the first case we retain the 5th 
order contribution to the mass matrix, while in the second we have a unique choice for 
the 2 nd family and the higgs, i.e., H = F 2R and F 2R accommodates the right-handed 
fields of the second generation. 

In order to calculate the contribution of the operator in eq. (p8|) to the mass matrices, 
we should properly contract the various fields involved in the NR-term. In principle, 
the numerical coefficient in front of the desired operator is a linear combination of the 
Clebsch Gordan coefficients presented in the Table, each of them multiplied by a dif- 
ferent phase factor. Our ignorance about the numerical coefficients of the mass matrix 
entries has been minimized in the unknown phase factors. However, we can make a 
natural assumption that the biggest contributions come from contractions occurring 
first for the fields belonging to the same sector. Recalling now that F 2R and H = F^r 
originate from the same sector b 2 , while F^l, H = F^ R and h\ 2 are obtained from b^. 
We find that 

(F 4L H)(F 2R H)h 12 - G 
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i.e., this operator contributes to down quark and charged lepton mass matrices. Since 
this contribution is the second largest after the tree level term F iL F 5R h 12 , we identify 
(|99|) with the 23 - entry of the corresponding mass matrices. It is clear therefore, that in 
this picture F^r, F 2 r accommodate the right components of the 3 rd and 2 nd generations 
respectively, whereas F^ contains the left fermions of the heavy generation. In order 
to obtain non-zero s-quark and fi- masses, we need to fill in the 32 entry of the down 
and charged lepton mass matrices with a higher order operator, so that the 2x2 lower 
block of the corresponding mass matrices exhibits a structure of the asymmetric type 
considered in the previous section, 



where rjij stand for higher order NR- contributions, while Fn represents in general one 
of the two remaining left-handed fourplets -F(i.3)l- To determine which of the latter will 
accommodate the 2 nd generation and calculate 77's, one has to proceed above the fifth 
order and find the relevant non-vanishing correlators. For example, choosing a new flat 
direction in which £ 2 , $2 singlets develop non-zero VEVs while interpreting F 3 l, F 3 r as 
the third generation, ^ 22 may arise from the fifth order NR-operator -Fil-Fir^G^ i n 
Eq. (|96|) . Interestingly, this is a n = operator according to our classification, however it 
is suppressed compared to a tree level term due to the presence of the 'effective' flavour 
factor 5 2 = f5l $2> . Furthermore, higher NR-terms will certainly involve n = 2, 3 etc 



operators. Thus, it is clear that the above procedure will require n > 2 operators and 
the analysis will be more involved than the field theory model described in the our 
earlier sections. 

As a matter of fact, a detailed analysis requires also the examination of all possible 
flat directions as well as calculation of the non-renormalizable contributions to even 
higher orders, since one has to ensure that the necessary Weinberg-Salam doublets 
living in some combination of our /z 3 , /i 12 bidoublets remain massless up to this order. 
For the moment, in our first approach to this model, we have been able to show that 
the rather complicated string construction stays in close analogy with the field theory 
approach presented in the early sections. 




(100) 



string 
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10 Conclusions 



We have examined Yukawa textures within a string inspired SU(4) x 0(4) model 
extended by a gauged U(l)x family symmetry and non-renormalisable operators above 
the unification scale of the form in Eq.|27|. These operators factorise into a factor {HH) 
and a factor involving the singlet fields 6, 9. The singlet fields 8, 9 break the U(l)x 
symmetry and provide the horizontal family hierarchies while the H, H fields break the 
SU(4)®SU(2)l®SU(2)k symmetry and give the vertical splittings arising from group 
theoretic Clebsch relations between different charge sectors. The factor {HH) also 
provides an additional flavour independent suppression factor S which helps the fit. 
The quark and lepton masses and quark mixing angles are thus described at high 
energies by a single unified Yukawa matrix whose flavour structure is controlled by a 
broken U(l)x family symmetry, and all vertical splittings controlled Clebsch factors. 
An important feature of the scheme is the existence of Clebsch zeroes which allow an 
entirely new class of textures to be obtained. For example the RRR solutions 3 and 
5 may be reproduced by this scheme which are complementary to the RRR solution 2 
favoured by the IR approach. 

In addition to the symmetric textures we have also performed a completely new 
analysis of the non-symmetric textures which are motivated by the string construction. 
A global fit to the fermion mass spectrum with 3 dof is described, in which three models 
in Table 3 are singled out with x 2 /dof< 1.Q At this level of difference of y 2 between 
models, the \ 2 test is subject to large statistical fluctuations. Therefore, we do not 
statistically distinguish between the fits in Tables 3,§ since both contain good fits to 
the data with x 2 /dof< 1. However we have a theoretical preference for the models 
in Table 3 since these models result from the operators in Eq.^7| where the family 
hierarchies are accounted for by the U(l)x symmetry, as explained in section |7|. By 
contrast, the models in Table || result from the operators in Eq.[26| and are essentially 
an up-dated version of those previously considered in ref . |18| . 

The string analysis performed in the later sections of the paper lends some support 

to the approach followed in this model. In the string model, the U(l) family symmetries 

14 By comparison a recent paper |^9) performed a global \ 2 analysis for some 50(10) models, 
including the mass and mixing data. With 3 dof, they obtain a x 2 /dof~ 1/3 for the best model. 
While our fit to model 3 in Table 3, for example, has a smaller x 2 /dof than this, it is difficult to make 
a comparison as in ref. |^9| quark mass correlations from data, as well as the effect of large tan (3 on 
mi, has been included. Also note that these involve the soft terms, thus a larger number of parameters 
are involved in the fit. 
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are a consequence of the string construction, but there are four of them with one being 
anomalous. There are several singlets (charged under the family group) to take the role 
of the 9 fields and the n = 1 operators involving a factor of HH are clearly expected 
in the effective theory below the string scale. We have shown that operators such as 
G which were simply pulled out of thin air in the earlier parts of the paper may in 
fact originate from string theory. As an example we constructed explicitly the lower 
2x2 block in Eq. |100| which has the characteristic asymmetric structure of the Yukawa 



textures considered earlier. It will be noted however that the lower 2x2 block in 
Eq. |100| does not correspond precisely to the ansatz in Eq.|35|. Within the given string 



construction, such an ansatz does not appear to be possible. The reason is the extra 
U(l) symmetries and the other discrete-like symmetries (selection rules) left over in the 
low energy model. A new string construction with a new boundary condition on the 
string basis is required in order to make contact with the phenomenologically preferred 
ansatze. This will be the subject of future work. 
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Appendix 1. n = 1 Operators 



The n = 1 operators are by definition all of those operators which can be constructed 
from the five fields FFhHH by contracting the group indices in all possible ways, as 
discussed in Appendix 1. After the Higgs fields H and H develop VEVs at M GUT of the 
form (H ab ) = (H A1 ) = uh, (H ax ) = (Ha) = the operators listed in the appendix 
yield effective low energy Yukawa couplings with small coefficients of order Mq UT /M 2 . 
However, as in the simple example discussed previously, there will be precise Clebsch 
relations between the coefficients of the various quark and lepton component fields. 
These Clebsch relations are summarised in Table [7|, where relative normalisation factor 
has been applied to each. The table identifies which SU(4) and SU(2) structures have 
been used to construct each individual operator by reference to Eqs. |103| . |104 . 



The n = 1 operators are formed from different group theoretical contractions of the 
indices in 

O og* = F aa Fp x hlH JZ H pw . (101) 

It is useful to define some SU(4) invariant tensors C, and SU(2)# invariant tensors R 
as follows: 

(d)? = $ 
(c 6 )5 = 

(C 10 )^ = sffi + sffl 

(RlTy = ^ 

{R,r y : = 6*8:- ~8*5™, (m 

where 8p, e a /3w X , 8®, e wz are the usual invariant tensors of SU(4), SU(2)^. The SU(4) 
indices on Ci 5 6,io,i5 are contracted with the SU(4) indices on two fields to combine them 
into 1, 6, 10, 15 representations of SU(4) respectively. Similarly, the SU(2)r indices on 
i?i.3 are contracted with SU(2)r indices on two of the fields to combine them into 1, 3 
representation of SU(2)#. 

The SU(4) structures in Table |^ are 

i- (Ci)£(d)z 

II. (C 15 )^(C 15 )^ 
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Table 7: When the Higgs fields develop their VEVs, the n = 1 operators lead to the effective 
Yukawa couplings with Clebsch coefficients as shown. 
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and the SU(2) structures are 



III. 
IV. 
V. 
VI. 



I. 
II. 
III. 

IV. 
V. 
VI. 



(c B )£(c 6 )£ 

(Gio)^(C 10 )S 
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(Ci6)S(Cib)£, 



3 J«ir 
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e xz e 
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(Rx)y(Ri 



XT 

wq' 



(103) 



(104) 



The operators are then given explicitly by contracting Eq. |101 with the invariant tensors 
of EqJlQ2| given by Table and Eqs. |103|J104 . 
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Appendix 2. Review of Analysis of Ref . [|18 



In ref. [18 we assumed that the Yukawa matrices at Mx are all of the form 



/ 0(e 2 ) 0(e 2 ) 

X U,D,E,N = 0(e 2) 0(e) 0(e) I ^ ( 105) 

V 0(e) 0(1) 

where e << 1 and some of the elements may have approximate or exact texture zeroes 
in them. First, we examine closer the assumption that the operator in the (33) position 
of the Yukawa matrices is the renormalisable one. It has been suggested in the past 
that the large value of tan j3 required by the constraint 

\{M GU t) = h{M GUT ) = X T {M GUT ) (106) 

such as is predicted by the renormalisable operator, leads to some phenomenological 
problems. One such problem is that a moderate fine tuning mechanism is required 
to radiatively break the electro-weak symmetry in order to produce the necessary 
hierarchy of Higgs VEVs t>i/t>2 ~ m t /mb []3~0"1 , [j31jl . One could set about trying to 



extend the present model in a manner that would lead to an arbitrary choice of tan /3, 
for example by introducing extra Higgs bidoublets. This route has its disadvantages in 
that a low value of tan f3 has been shown |32[ in most schemes to be inconsistent with 
Xb^Mcur) = K(Mcut) unification if the tau neutrino mass constitutes the hot dark 
matter requiring the Majorana mass of the right handed tau neutrino to be M^ T ~ 
O(10 12 ) GeV. To a very good approximation, the largest diagonalised Yukawa coupling 
in X 1 is equal to its 33 entry A33. (One may obtain small tan/? solutions consistent 
with mb-m T unification and an intermediate neutrino scale, in specific models: Either 



large mixing in the /i — r charged leptonic sector has to occur [ 33 1 or the Dirac-type 



Yukawa coupling of the neutrino has to be very suppressed fl34|.) 

To force things to work in a generic scheme, one solution could be to use a non- 
renormalisable operator in the 33 position which has some Clebsch factor x > 1 such 
that 

X t (M GUT ) = xX b {M GUT ) = xX T (M GUT ). (107) 



Eq. |107| would preserve the bottom-tau Yukawa unification, but lower the prediction 
of tan/3 due to the bigger contribution to the top Yukawa coupling. It may only 
be reasonable to examine n = 1 operators in this context since we know that the 
third family [H| Yukawa coupling is ~ 0(1) and higher dimension operators could 
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be expected to provide a big suppression factor. Systematically examining the n = 1 
operators we find that only the operator 33 , which leads to the prediction 



A t (M ( 



GUT) 



2X b (M ( 



GUT) 



2A T (M« 



GUT) 



(108) 

can decrease tan (3. The change is minimal, from 56.35 to 55.19 for ois(Mz) = 0.117 and 
M^ T = O(10 12 ) GeV. The reason that the change is minimal is due to the fact that the 
Yukawa couplings are approximately at their quasi fixed points [[35] and so even a large 
change to Xt,b,r(Mx) produces only a small change in A 4j b iT (m i ), which are the quantities 
that require a high tan/? through the relations in Eq.^T]. Another possibility would 
be to include 0^,0 33 which would allow arbitrary tan/5 (in particular intermediate 
tan/3 ~ 10 — 20.) However, this would reduce the predictivity of the scheme as tan/? 
would become an input. One might also be skeptical about whether a parameter 
~ 1 could be generated by a non-renormalisable operator in a perturbative scheme. It 
would certainly require the heavy mass scales M to be very close to the VEVs H, H, 8, 9 
and we might therefore naively expect large corrections to any calculation based on 
this model. We thus abandon these ideas and continue with the usual renormalisable 



operator in the 33 position of the Yukawa matrices that leads to Eq. |106| . We note in 
any case that a recent analysis explains that in gauge mediated supersymmetry 
breaking models, the radiative mechanism of electroweak symmetry breaking can be 
such that no fine tuning occurs for large tan /3. In these models high tan (3 admits 
solutions of the hot dark matter problem in which the Yukawa couplings unify ||32|| . 

The hierarchy assumed in Eq. |105| allows us to consider the lower 2 by 2 block of the 
Yukawa matrices first. In diagonalising the lower 2 by 2 block separately, we introduce 
corrections of order e 2 and so the procedure is consistent to first order in e. We found 
several maximally predictive ansatze that were constructed out of the operators whose 
Clebsch coefficients are listed in table |4] for the n = 1 operators. The explicit n = 1 
operators in component form are listed in the Appendix 1. We label the successful 
lower 2 by 2 ansatze Af 



A l = 


og-og 
og o 33 


(109) 


A 2 = 


o og-og' 
_ og 33 


(110) 


A 3 = 


o og-og' 
_ og 33 


(111) 


A A = 


o og ' 

Og - Og O33 


(112) 
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O33 

nG 

^23 
O33 



o£ o 



^23 
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(113) 
(114) 
(115) 
(116) 



We now note that solutions A 2 -s require a parameter if 23 ~ 0(1) to attain the correct 
and V c b- Any calculation based on the hierarchy assumed in Eq. |105| is therefore 
inconsistent and so we discard these solutions. We also note that O32 only has the 
effect of fixing V c b to a good approximation and so can consist of any operator in 
Table [7] that has a different Clebsch coefficient for up quark and down quark Yukawa 
couplings. The precise operator responsible for V c b has no bearing on the rest of the 
calculation and we therefore just make an arbitrary choice of Og for the rest of this 
paper. We also note that for the phenomenologically desirable and predictive relation 

Xg(M GUT ) 



GUT) 



(117) 



to hold, we may replace Og - Og in A x with 0% 2 + Og, Og + Og or any other 



combination of two operators which preserves Eq. |117| and allows Xg to be smaller and 

D E 

independent of A 2 2 • In fact, the preferred solution is that the dominant operator in 
that position be Og which does not give a contribution to the up quark mass. Then, 
a subdominant operator would be responsible for the entry Xg and would therefore be 
suppressed naturally by one or more powers of e. 



Appendix 3. Flatness Conditions in the String Model. 



We give here the constraints on the various singlet VEVs obtained from the F 
and D flatness conditions in the string spectrum of the model in section 8. From the 
F-flatness of the superpotential one derives 18 conditions, which are the following : 



66 + C1 2 + cl 
66 + C? + C 2 2 










lis) 
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tei + C1C1 + C2C2 

C1C2 + C1C2 

2^T2Ci + 2 $ sCi + $4(2 
2<%Ci + 2 $ 3Ci + $4(2 
2$ n < 2 + 2*sC 2 + $ 4Cl 

2$T2C2 + ^$3C2 + $4Cl 

$126 + ^$36 

$12^4 + ^$ 3 fl 
$12& + 2^s6 
$ 12& + ^$3& 
$ 12& + ^36 
^126 + 
$126 + ^$3& 
$126 + ^$3& 





















(119) 



Now, a possible choice of non- zero singlet VEVs which satisfy the system ( |119| ), is 

< $ 12 > , < i>^ 2 > , < a > , < 6 > ^ (120) 

accompanied with non - zero VEVs of the following two hidden (octets under SU(8)h) 
fields 

< Z 5 > , < Z! 3 > ^ (121) 
Taking all other singlet and hidden field VEVs equal to zero, the D- flatness conditions 



read |]39| 



2|$ 



12 1 



N 2 -I6l 2 + 



3a M 
2tt 



2^ 5 | 2 -I6l 

2ie 1 r-i6i 2 -2|«i>r 2 r + ^3l 2 + l^i 

2|$ 12 | 2 + |£ 1 | 2 -^|Z 



/ 12 










(122) 
(123) 
(124) 
(125) 
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The scale of the non- zero singlet VEVs is determined by the above conditions. There 
are five equations to determine four parameters, thus one has the freedom to fix one of 
the non zero VEVs in Eqs. (|120|J121|) from phenomenological requirements. In any case, 
from the above equations it turns out that the natural scale of the non zero VEVs are 
of the order ^-Mpi. For a u ~ 10 _1 one can see that their magnitude is of the required 
order to contribute in the mass operators. 



50 



References 

[1] M. S. Chanowitz, J. Ellis and M. K. Gaillard, Nucl. Phys. B128 (1977) 506; A. J. 
Buras, J. Ellis, M. K. Gaillard and D. V. Nanopoulos, Nucl. Phys. B135 (1978) 
66; M. B. Einhorn and D. R. T. Jones, Nucl. Phys. B196 (1982) 475; J. Ellis, D. 
V. Nanopoulos and S. Rudaz, Nucl. Phys. B202 (1982) 43. 

[2] B. C. Allanach and S. F. King, Phys. Lett. B353 (1995) 477. 

[3] B. Ananthanarayan, G. Lazarides and Q. Shan, Phys. Rev. D44 (1991) 1613; H. 
Arason et al, Phys. Rev. Lett. 67 (1991) 2933; L. Hall, R. Rattazzi and U. Sarid, 
Phys. Rev. D50 (1994) 7048. 

[4] H. Fritzsch, Phys. Lett. B70 (1977)436. 

[5] H. Georgi and C. Jarlskog, Phys. Lett. B86 (1979) 297; S. Dimopoulos, L. Hall 
and S. Raby, Phys. Rev. D45 (1992) 4192. 

[6] P. Ramond, R. Roberts and G. G. Ross, Nucl. Phys. B406 (1993) 19. 

[7] J. Harvey, P. Ramond and D. Reiss, Nucl. Phys. B199 (1982) 223; Phys. Lett. 
B92 (1980) 309; C. Wetterich, Nucl. Phys. B261 (1985) 461; Nucl. Phys. B279 
(1987) 711; J. Bijnens and C. Wetterich, Phys. Lett. B176 (1986) 431; Nucl. 
Phys. B283 (1987) 237; Phys. Lett. B199 (1987) 525; F.J. Gilman and Y. Nir, 
Ann. Rev. Nucl. Part. Sci 40 (1990) 213; S. Dimopoulos, L. J. Hall and S. Raby, 
Phys. Rev. D45 (1992) 4192; Phys. Rev. Lett. 68 (1992) 1984; Y. Achiman and 
T. Greiner, Nucl. Phys. B443 (1995) 3; E. Papageorgiu, Z. Phys. C64 (1994) 509; 
Z. Phys. C65 (1995) 135; Phys. Lett. B343 (1995) 263; Y. Grossman and Y. Nir, 
Nucl. Phys. B448 (1995) 30; H. Fritzsch and Zhi-zhong Xing, Phys. Lett. 353 
(1995) 114; C. H. Albright and S. Nandi, Mod. Phys. Lett. All (1996) 737; Phys. 
Rev. D53 (1996) 2699. 

[8] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B147 (1979) 277. 

[9] L. Ibanez and G. G. Ross, Phys. Lett. B332 (1994) 100. 

[10] H. Dreiner, G. K. Leontaris, S. Lola, G. G. Ross and C. Scheich, Nucl. Phys. B436 
(1995) 461. 



51 



[11] G. G. Ross, Phys.Lett. B364 (1995) 216; B. C. Allanach and S. F. King, RAL 
TR 97-007, SHEP 97-01, frep-ph/9703295 . 

[12] G. Anderson, S. Raby, S. Dimopoulos, L. Hall and G. Starkman, Phys. Rev. D49 
(1994) 3660. 

[13] R. Barbieri and L. Hall, Phys. Lett. B377 (1996) 76; R. Barbieri and L. J. Hall, 
LBL-38381, [hep-ph/9605224] 



[14] Z. Kakushadze and S. H. Henry Tye, CLNS-96-1433, |hep-th/9610T06 

[15] J. Pati and A. Salam, Phys. Rev. D10 (1974) 275. 

[16] I. Antoniadis and G. K. Leontaris, Phys. Lett. B216 (1989) 333. 

[17] I. Antoniadis, G. K. Leontaris and J. Rizos, Phys. Lett. B245 (1990) 161; G.K. 
Leontaris, Phys. Lett. B372(1996)272. 

[18] S. F. King Phys. Lett. B325 (1994) 129; B. C. Allanach and S. F. King, 
Nucl.Phys.B456(1995)57; ibid B459(1996)75. 

[19] I. Antoniadis, C. Bachas, and C. Kounnas, Nucl.Phys.B289 (1987) 87; I. An- 
toniadis and C. Bachas, Nucl.Phys.B298 (1988) 586; B. Greene, K.H. Kirklin, 
P.J. Miron and G.G. Ross, Nucl.Phys.B292 (1987) 606; A. Font, L.E. Ibanez, F. 
Quevedo and A. Sierra, Nucl.Phys. B331 (1990) 421; I. Antoniadis, J. Ellis, J. 
Hagelin, and D.V. Nanopoulos, Phys. Lett. B231 (1989) 65; A. Faraggi, Phys. 
Lett. B278 (1992) 131; Nucl.Phys.B387 (1992) 239. 

[20] B. Greene, K.H. Kirklin, P.J. Miron and G.G. Ross , Nucl.Phys.B292 (1987) 606; 
A. Font, L.E. Ibanez, F. Quevedo and A. Sierra, Nucl.Phys. B331 (1990) 421; 
D. Bailin, A. Love and S. Thomas, Nucl.Phys.B298 (1988) 75; J.A. Casas, E.K. 
Katehou and C. Munoz, Nucl.Phys.B317 (1989) 171. 

[21] I. Antoniadis, J. Ellis, J. Hagelin, and D.V. Nanopoulos, Phys. Lett. B231 (1989) 
65; A. Faraggi, D.V. Nanopoulos and K. Yuan, Nucl.Phys.B335 (1990) 347. A. 
Kagan and S. Samuel, Phys. Lett. B284 (1992) 289; J. Lopez, D.V. Nanopoulos, 
and K. Yuan, Nucl.Phys.B399 (1993) 654; A. Faraggi, Phys. Lett. B278 (1992) 
131; Nucl.Phys.B387 (1992) 239. S. Chaudhuri, G. Hockney, and J. Lykken, [Kep^l 
th/9510241] ; G.K. Leontaris, Phys. Lett. B372(1996)272, |hep-ph/9601337| . 



[22] A. E. Farragi and J. C. Pati, |hep-ph/ 9703235 . 



52 



[23] G. K. Leontaris, S. Lola, C. Scheich and J. D. Vergados, Phys. Rev. D53 (1996) 
6381. 



[24] U. Amaldi et al, Phys. Rev. D36 (1987) 1385; P. Langacker and M.Luo, Phys. 
Rev. D44 (1991) 817; J. Ellis,S. Kelley and D.V. Nanopoulos, Phys. Lett. B249 
(1990) 441; Nucl. Phys. B373 (1992) 55; U.Amaldi, W. de Boer and H.Fustenau, 
Phys. Lett. B260 (1991) 447; C. Giunti, C.W. Kim and U.W. Lee, Mod. Phys. 
Lett. A6 (1991) 1745; H. Arason et al., Phys. Rev. D46 (1992) 3945; F. Anselmo, 
L.Cifarelli, A.Peterman and A.Zichichi, Nuovo Cimento 105A (1992) 1179; P. 
Langacker and N. Polonski, Phys. Rev. D47 (1993) 4028; A.E. Faraggi and B. 
Grinstein, Nucl. Phys. B422 (1994) 3. 

[25] Particle Data Book, Phys. Rev. D54, (1996) 1. 

[26] L. J. Hall, R. Rattazzi and U. Sarid, Phys. Rev. 50D (1994) 7048; M. Carena, M. 
Olechowski, S. Pokorski and C. E. M. Wagner, Nucl. Phys. 426B (1994) 269. 

[27] P. Binetruy and P. Ramond, Phys.Lett. B350 (1995) 49; P. Binetruy, S. Lavignac 
and P. Ramond, |hep-ph/9601243| ; P. Ramond, [hep-ph/9604251| , Workshop on 



Frontiers in Quantum Field Theory, Osaka, 1995; S. Lavignac, [hep-ph/ 9605355 



Proceedings of 31st Rencontres de Moriond, Les Arcs, 1996; S. Lavignac, 



hep-ph/96 10257| , NATO Advanced Study Institute on Masses of Fundamental 



Particles, Cargese, France, 1996. 

[28] M. Dine, N. Seiberg and E. Witten, Nucl.Phys.B289 (1987) 585. 

[29] T. Blazek, M. Carena, S. Raby and C. Wagner, OHSTPY-HEP-T-96-014, 
[hep-ph/9608273 . 



[30] R. Rattazzi and U. Sarid, Phys. Rev. D53 (1996) 1553. 

[31] M. Carena, M. Olechowski, S. Pokorski and C.E.M. Wagner, Nucl. Phys. B426 
(1994) 269. 

[32] B.C. Allanach and S.F. King, Phys. Lett. B353 (1995) 477; F. Vissani and A.Y. 
Smirnov, Phys. Lett. B341 (1994) 173; A. Brignole, H.Murayama and R.Rattazzi, 
Phys. Lett. B335 (1994) 345. L.J.Hall, R.Rattazzi and U.Sarid, Phys. Rev. D50 
(1994) 7048. 

[33] G. K. Leontaris, S. Lola and G. G. Ross, Nucl. Phys. B 454 (1995) 25. 



53 



[34] S. Dimopoulos and A. Pomarol, Phys. Lett. B 353 (1995) 222. 
[35] B.Schrempp Phys.Lett. B344 (1995) 193. 

[36] M. Dine, Y. Nir and Y.Shirman, SCIPP 96/30, WIS-96/29, WIS-96/29/Jul-PH, 
[hep-pri/9607397 . 

[37] A. Font, L. Ibanez H.P. Nilles and F. Quevedo, Phys. Lett.B210 (1988)101; 
L. Dixon, D Friedan, E. Martinec and S. Shenker, Nucl. Phys. B282(1987)13; 
J. Rizos and K. Tamvakis, Phys. Lett. B262 (1991) 227; 
S. Kalara, J. Lopez and D.V. Nanopoulos, Nucl.Phys.B 353 (1991) 650 
A. E. Faraggi, Nucl. Phys. B 403(1993) 101. 

[38] I. Antoniadis, G.K. Leontaris and N.D. Tracas, Phys. Lett.B279(1992)58; 
K. R. Dienes, A. Faraggi and J. March-Russell, Nucl. Phys.B476(1996)44. 

[39] G.K. Leontaris and N.D. Tracas, Zeit. Phys. 56 (1992) 479. 



54 



